r 

ARMY  ENSINCEH  MTEMAVS  CXPtRIHCNT  STATION  VICKSBUM~£TC  r/s  , 

1  "*vt  ""'""'ox  cogc'STiRnffM^&u. 

UNCLASSIFIED  WES/MP/SL-A1-25  ^ 

1 

m 

END 

10*81 

DT1C 

■ 

4  i 

■■■ 

BB 

p w 

1  -  '  ~  1 

mm 


MISCELLANEOUS  PAPER  SL-8L25 


A  ONE-DIMENSIONAL  PLANE  WAVE 
PROPAGATION  CODE  FOR  LAYERED 
RATE-DEPENDENT  HYSTERETIC  MATERIALS 


John  O.  Curtis 
Structural  Laboratory 

U.  S.  Army  Enginaar  Waterways  Experiment  Station  f  ’’ 

P.  O.  Box  631,  Vicksburg,  Miss.  39180  LJ  *  1^“* 

^ELECTE 

September  1981  V  ^  OCT  8  198t 

Final  Raport 

I  Apvrosxl  For  Public  RUmm;  Distribution  UnUmltod  I 


<##Sbw»«ir ; 


iw  Assistant  Secretary  of  the  Army  (R6>D) 
Department  of  the  Army 
Washington,  D.  C.  20315 

(Me  Project  No.  4AI6II0IA9ID 


41  10  7  08» 


to 

to 


Unclassif led _ 

security  CLASSIFICATION  of  this 


AGE  Datm  Kntered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  _ J  2.  GOVT  ACCESSION  NO. 

Miscellaneous  Papeip  SL-81-25  /  ft  ^ 

3.  RECIPIENT'S  CATALOG  NUMBER 

mvjzL 

4  TITLE  (mnd  Subtitle) 

A  QNE-MMENS IONAL  PLANE  WAVE  PROPAGATION  CODE  FOR 

Jayered  ^te-^ependent^Jysteretic  materials* 

b^?Vpe  of  report  A  PERIOD  covered 

Final  repesrt  _ ...  - 

t  PERFORMING  ORG.  REPORT  NUMBER  j 

n  ' '  7  ^  ~  ^  i 

7.  AUTHORS 

John  0. /Curtis 

B.  ’CONTRACT  OR  GRANT  NUUBEH(m)  J 

S-  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

U.  S.  Army  Engineer  Waterways  Experiment  Station 
Structures  Laboratory 

P.  0.  Box  631,  Vicksburg,  Miss.  39180 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Pro jeAt^  4A1611|4lA91D^ 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  /' 

Assistant  Secretary  of  the  Army  (R6D)  /  ' 

Department  of  the  Army 

Washington,  D.  C.  20315 

U.  WBABT  DATE 

September  J.981/ 

13.  NUMBER  OF  PAGES 

90 

<4.  MONITORING  AGENCY  NAME  «  ADDRESS?/?  dlHoront  tram  Controlling  Ollier) 

- V.  / 

:■/  /'  ■  J 

is.  SECURITY  CLASS,  (ol  Olio  rmport) 

Unclassified 

ls«.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 

l«.  DISTRIBUTION  STATEMENT  (ol  Olio  Rmport) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tho  that  root  ontorod  In  Block  20,  It  dlltoront  Irom  Report) 

i».  supplementary  notes 

Available  from  National  Technical  Information  Service,  5285  Port  Royal  Road, 
Springfield,  Va.  22151. 

If-  KEY  WORDS  (Continue  on  f*rtrM  aid*  It  nacMMrr  and  identify  by  block  number) 

Compaction  (Soils)  Hysteretic  media 

Constitutive  models  0NED3P  (Computer  program) 

Finite  element  method  Wave  propagation 

20.  ABSTRACT  fCarftet  an  rarer  mm  ft  neemmeery  mad  Identify  by  block  number) 

^  Recent  experimental  data  on  soils  showing  that  the  application  of  loads 
with  submillisecond  rise  times  results  in  significant  rate-dependent  com¬ 
pressibility  behavior  have  prompted  the  need  for  a  one-dimensional  plane  wave 
propagation  code  for  layered  rate-dependent  hysteretic  or  visco-compacting 
materials.  -Accordingly;?  gn  explicit  one-dimensional  finite  element  code- 
named  0NED3^  h«s--be*i^devdioped  which  incorporates  a  three-parameter  (spring 

_  -  io  (Continued)  ^/_ 

00  ,'STn  WZ  I'O'-na.o,  ■  no  v  obsolete  Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  FACE  (When  Dmtm  Entered) 


Unclassified  _ _ 

ItCUHITY  CIAMIHCAT1QN  Or  THU  PAOt(TWmi  *>«<«  tnlmnd) 

ABSTRACT  (Continued 1 . 

and  dashpot)  mechanical  model  consisting  of  a  linear  spring  and  dashpot  in 
series  coupled  in  parallel  with  a  piecewise  linear  hysteretic  or  compacting 
spring. 

0NED3P  allows  for  th&  analysis  of  a  column  of  multilayered  viscocom¬ 
pacting  soils  loaded  by  a  digitized  surface  pressure-time  history.  Any  set  of 
consistent  units  may  be  used  with  this  code  and  results  may  be  obtained  in  the 
form  of  stress-,  strain-,  acceleration-,  velocity-,  and/or  displacement-time 
histories  as  well  as  stress-strain  curves  using  standard  Calcomp  software. 

Several  demonstration  problems  were  calculated  using  0NED3P  to  evaluate 
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A  ONE-DIMENSIONAL  PLANE  WAVE  PROPAGATION 
CODE  FOR  LAYERED  RATE-DEPENDENT 
HYSTERETIC  MATERIALS 


PART  I:  INTRODUCTION 


Background 

1.  Material  property  test  data  recently  acquired  in  the  Geo¬ 
mechanics  Division  of  the  Structures  Laboratory  at  the  U.  S.  Army  Engi¬ 
neer  Waterways  Experiment  Station  (WES)  have  revealed  that  the  dynamic 
compressibility  response  of  various  soils  subjected  to  loadings  with 
submillisecond  rise  times  is  both  qualitatively  and  quantitatively  dif¬ 
ferent  than  their  response  to  slower  loadings.^  Reference  1  cites  sev¬ 
eral  examples  of  test  results  in  which  the  stress-strain  response  of 
soils  was  much  stiffer  during  rapid  loading  conditions  than  for  quasi¬ 
static  experiments. 

2.  Furthermore,  field  measurements  recently  acquired  during 
shallow-buried  structures  experiments  (References  2  and  3)  indicated 
that  for  surface  loading  rise  times  on  the  order  of  0.01  to  0.1  msec, 
high-amplitude  (10  to  40  MPa)  stress  waves  traveled  faster  through  the 
sand  cover  above  the  structures  than  would  be  predicted  from  seismic  ve¬ 
locity  data  or  from  uniaxial  strain  test  data  generated  in  the  laboratory 
using  loadings  with  rise  times  on  the  order  of  a  few  milliseconds. 

Purpose  and  Scope 

3.  The  purpose  of  this  report  is  to  describe  the  development  and 
evaluation  of  a  one-dimensional  plane  stress  wave  propagation  code 
called  0NED3P  which  treats  layered,  nonlinear,  rate-dependent,  hysteretic 
materials . 

4.  A  description  of  the  constitutive  relationship  used  in  0NED3P, 
which  is  represented  by  a  three-parameter  mechanical  model,  is  given  in 


3 


Part  II.  Part  III  contains  a  description  of  the  major  features  of 
0NED3P,  including  its  solution  algorithm  and  its  treatment  of  different 
boundary  conditions.  The  capabilities  of  the  code  are  checked  against 
available  analytical  solutions  and  other  code  calculations  in  Part  IV. 
Part  V  describes  how  the  model  parameters  can  be  evaluated  from  labora¬ 
tory  and  field  data  and  goes  on  to  compare  0NED3P  calculation  results 
with  those  from  field  experiments.  Finally,  a  user's  guide  is  presented 
in  Appendix  A. 
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V. 


PART  II:  MODEL  DESCRIPTION 


Mechanical  Model 


5.  In  general,  a  rate-dependent  constitutive  model  should  relate 
stress,  strain,  stress  rate,  and  strain  rate  in  the  following  functional 
form: 


a  =  f(o,  e,  e)  (1) 

where  the  dot  indicates  time  differentiation. 

6.  The  literature  (References  A  and  5)  shows  that  numerous  rela¬ 
tionships  among  forces,  displacements,  loading  rates,  and  velocities  may 
be  written  by  developing  the  governing  equations  for  various  combina¬ 
tions  of  linear  mechanical  elements;  namely,  springs  and  dashpots.  The 
springs  generate  forces  (or  stresses)  proportional  to  displacements  (or 
strains) ,  while  the  dashpots  generate  forces  proportional  to  velocity 

(or  strain  rate).  An  example  of  a  mechanical  model  whose  governing  equa¬ 
tion  looks  like  Equation  1  but  which  is  still  relatively  simple  to  work 
with  is  shown  in  Figure  la.  Furthermore,  Reference  1  has  already  demon¬ 
strated  that  such  a  model  can  be  used  to  simulate  the  rate-dependent 
stress-strain  reponse  of  a  single  particle. 

7.  -It  can  be  shown  that  the  equation  which  governs  the  behavior 

4 

of  this  three-parameter  model  is 

M,  MM 

a  +  —  a  =  - -  e  +  (M„  +  M.  )e  (2) 

n  n  2  1' 

The  M^  and  M2  functions  in  Equation  1  bound  the  stress-strain  behav¬ 
ior  of  the  material  in  the  following  ways.  First,  for  very  slow  loading 
rates  (0  ,  e  -*■  0) ,  M^  describes  the  complete  stress-strain  response 
and,  hence,  represents  the  "static"  behavior  of  the  material.  On  the 
other  hand,  for  extremely  fast  loading  rates,  the  dashpot  acts  as  though 
it  were  rigid  and  the  sum  of  M^  and  M2  governs  the  model  behavior. 


j  <rlt) 


Mx  =CONSTANT 
ff  =CONSTANT 


a.  ARRANGEMENT  OF  SPRINGS  AND  DASHPOT 
IN  MECHANICAL  MODEL 


£ 


b.  STATIC  STRESS-STRAIN  CURVE  DEFINING  THE  M,  FUNCTION 

Figure  1.  Details  of  the  proposed  visco¬ 
compacting  mechanical  model  for  soils 
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Therefore  +  M2  represents  the  upper  bound  of  material  stiffness. 

Any  material  response  between  these  bounds  is  then  controlled  by  the 
value  assigned  to  n  . 

Finite  Difference  Form  of  the  Governing  Equation 

8.  Equation  1  is  written  for  constant  material  property  param¬ 
eters  ,  M2  ,  and  n  and  deals  with  total  values  of  stress  and 
strain  and  their  time  derivatives.  To  accommodate  nonlinear  properties. 
Equation  2  was  written  in  incremental  form  as 


Ao  +  Ao  =  M.  Ac  +  n  ( 1  +  7tMac 

M2  1  V  M2/ 


wherein  M^  ,  M2  ,  and  n  are  considered  to  be  constants  within  each 
increment  of  time.  In  fact,  as  a  first-order  approximation  to  obtaining 
a  model  of  the  visco-compacting  material  described  in  Part  I,  M2  and  n 
are  treated  as  constants  for  all  time.  M^  is  described  by  a  piecewise 
linear  stress-strain  curve  (with  an  arbitrary  number  of  segments)  which 
has  separate  loading  and  unloading  behavior  (Figure  lb) . 

9.  Using  the  subscript  i+1  to  refer  to  a  new  point  in  time,  i 
to  refer  to  the  present  time,  and  i-1  to  refer  to  the  previous  point 
in  time,  the  incremental  terms  in  Equation  3  may  be  written: 


Ao  =  o  .  , ,  -  o  . 
l+l  1 


Ac  -  ei+1  -  c± 

•  •  ai+l  ~  °i  °i  ~  °i-l 

Ao  °i+l/2  “  °i-l/2  At  At 

n  o 


£i+l  '  ei  Gi  ei-l 


where  At  =  t..,  -  t,  and  At  *  t.  t.  .  .  Substituting  Equation  4 
n  i+1  i  o  i  i-1 

into  Equation  3  yields  the  final  incremental  form  of  the 


constitutive  relationship  used  in  0NED3P: 


10.  As  long  as  remains  constant  over  any  two  consecutive 

time  increments,  AtQ  is  equal  to  At^  .  If,  however,  changes 

within  a  time  increment  the  following  technique  is  used  to 
evaluate  the  new  state  of  stress.  Consider  the  diagrams  shown 

in  Figure  2  where  the  strain  of  some  piece  of  the  material  has 
proceeded  from  point  a  to  point  b  to  point  d  on  two  consecutive  time  in¬ 
crements.  There  is  a  change  in  between  points  b  and  d.  0NED3P 

automatically  breaks  the  time  increment  from  b  to  d  into  two  (or  more, 
if  necessary)  subintervals  which  are  proportional  to  the  strain  segments 
(c^  -  and  (c^  -  e^) .  The  code  then  proceeds  as  shown  in  Figure  2  to 

calculate  the  new  stress  state  in  two  steps  (or  morn). 

11.  Other  techniques  for  handling  a  change  in  the  function 

were  tried  including  (a)  rewriting  Equation  3  to  include  variable 
values,  (b)  choosing  one  value  or  the  other  to  be  used  in  Equa¬ 
tion  5,  and  (c)  computing  an  average  value  of  over  an  interval. 

Assuming  that  there  were  no  programming  errors,  none  of  the  other  tech¬ 
niques  gave  results  as  consistently  stable  as  did  the  method  shown  in 
Figure  2. 


PART  III:  CODE  DESCRIPTION 


LJ.  Spatially,  GNED3P  employs  two-node  isoparametric  one¬ 
dimensional  finite  elements  in  which  displacement  is  assumed  ary 

linearly  between  the  nodes.  This  assumption  leads  to  a  constant  strain 
and,  hence,  constant  stress  within  each  element.  Furthermore,  half  of 
the  elements 's  mass  is  assigned  to  each  node  and  because  each  element 
is  assumed  to  Imvc  a  unit  cross-sectional  area,  the  stress  in  each  ele¬ 
ment  may  be  repiaicd  by  node  point  loads  (or  forces)  equal  in  magnitude 
to  the  element  stress.  Nodal  accelerations  are  found  by  summing 
the  forces  acting  on  each  node  and  dividing  by  the  node  mass. 

13.  This  scheme  allows  for  a  simple  visual  interpretation  of  what 
the  ON ED 3 P  code  denis  with.  Figure  3a  shows  a  column  of  continuous 
material,  the  horizontal  dimensions  of  which  have  no  meaning  (since  the 
code  works  with  a  unit  cross-sectional  area).  Figure  3b  shows  the 
equivalent  system  of  lumped  masses  and  mechanical  elements  with 

which  0NED3P  actually  works. 

14.  The  solution  algorithm  for  the  equivalent  system  is  shown 
graphically  in  Figure  4.  New  nodal  velocities  (V)  and  displacements  (d) 
are  calculated  by  a  simple  linear  integration  scheme: 

V  =  V  +  a  •  At  \ 

new  old  new  I 

>  (6) 

d  =  d  +  V  •  At  J 

new  old  new  / 


where  "a"  stands  for  acceleration.  A  mote  exact  integration  scheme 
was  tried  but  led  to  numerical  instabilities.* 

*  It  is  the  author's  belief  that  the  linear  integration  scheme  used  in 
0NED3P  serves  to  make  the  normally  lagging  response  of  an  element  (due 
to  the  finite  spacing  of  nodes  and  the  use  of  a  finite  time  increment) 
"catch  up"  to  the  true  solution  by  computing  larger  velocity  incre¬ 
ments  while  accelerating  and  smaller  increments  when  decelerating. 


i 

I 


f 
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* 


* 
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fl/G/0  Off  ffffff 

TRANSMITTING 

BOUNDARY 


a)  COLUMN  OF  CONTINUOUS  MATERIAL 


b)  EQUIVALENT  SYSTEM  FOR  ONED3P 


Figure  3.  0NED3P  interpretation  of  a  typical 
one-dimensional  problem 
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Figure  4.  Solution  algorithm  for  0NED3P 


Boundary  Conditions 


15.  Fixed  and  free  bottom  boundary  conditions  are  handled  in 
0NED3P  just  as  they  are  in  all  finite  element  codes;  namely,  the  acceler¬ 
ation  for  a  fixed-surface  node  is  always  set  equal  to  zero  while  a  free- 
surface  node  is  treated  like  any  other  node  within  the  material  column. 

16.  An  approximate  transmitting  bottom  boundary  has  also  been 
incorporated  into  0NED3P.  The  expression  for  particle  velocity  and 
strain  from  linear  elastic  one-dimensional  wave  propagation  theory  is:"* 
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V  =  Ce 


(7) 


where  C  is  the  wave  speed  which,  under  uniaxial  strain  conditions,  is 


C  = 


(8) 


M  being  the  constrained  modulus  and  p  the  density  of  the  material. 
To  implement  the  transmitting  boundary.  Equation  7  was  written  incre¬ 
mentally  as: 


V 

new 


V  ,  ,  +  C(e 
old  new 


Cold) 


(9) 


where  the  strain  at  the  boundary  node  was  taken  to  be  the  current  strain 
in  the  last  element. 

17.  The  wave  speed  is  computed  anew  at  each  time  step  as  a  func¬ 
tion  of  the  current  value  of  in  the  last  element  using  the  follow¬ 

ing  reasoning.  First,  the  tangent  modulus  at  any  point  on  the  dynamic 
stress-strain  curve  cannot  be  used  because  it  can  have  negative  values. 
Second,  the  current  value  of  was  found  to  be  insufficient  for  vis¬ 

cous  problems  in  which  the  dynamic  stress-strain  curve  was  much  different 
than  the  static  curve.  As  expected,  using  alone  resulted  in  a  soft 

boundary  for  highly  viscous  calculations.  Obviously  what  is  needed  is 
a  modulus  which  approaches  under  nonviscous  conditions  but  which 

takes  into  account  the  stiffer  dynamic  stress-strain  response  of  highly 
viscous  calculations.  One  measure  of  the  dynamic  response  in  a  material 
is  the  net  amount  of  energy  absorbed  at  a  given  point — in  other  words, 
the  area  under  the  stress-strain  curve.  Based  upon  these  observations, 
the  value  of  M  in  Equation  8  was  finally  taken  to  be: 

A, 

M  =  Mx  —  (10) 

s 

where 

A^  =  the  area  under  the  loading  portion  of  the  dynamic  stress- 
strain  curve 


I'kf'i  &1.  * 
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=  the  area  under  the  loading  portion  of  the  static  curve 
This  treatment  gives  stable  results  and  has  worked  quite  well  under  most 
conditions. 


Surface  Loading 

18.  Since  the  one-dimensional  column  being  simulated  by  0NED3P 
is  assumed  to  have  a  unit  cross-sectional  area,  force  and  stress  are 
synonymous  and  therefore  a  stress-time  history  may  be  applied  directly 
to  the  surface  node  as  a  force-time  history.  To  allow  for  generality 
of  input,  the  surface  forcing  function  used  in  0NED3P  must  be  digitized 
and  read  by  the  code  from  a  data  file.  With  only  a  minor  modification 
to  the  code,  a  velocity-time  history  could  be  applied  to  the  surface 
as  well. 


Consistent  Units 


19.  One  convenient  feature  of  0NED3P  is  that  any  set  of  consistent 
units  may  be  used  in  the  code.  Consistent  units  are  sets  of  units  which 
do  not  require  conversion  factors  to  make  calculations  balance  in  a  unit 
sense.  Any  of  the  sets  of  units  shown  in  Table  1  may  be  used  in  0NED3P 
to  generate  equivalent  results.  Note  that  set  D  represents  the  set  of 
units  normally  referred  to  as  SI  units. 

Plotting  of  Results 

20.  0NED3P  has  been  written  to  generate  several  types  of  plots 
using  standard  Calcomp  software.  Time  histories  of  stress,  strain, 
acceleration,  velocity,  and  displacement  may  be  obtained  as  well  as  plots 
of  total  dynamic  stress  versus  strain.  Further  information  on  plots  is 
available  in  Appendix  A. 
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Consistent  Sets  of  Units 


Energy  is  not  dealt  with  in  0NED3P ,  but  was  added  to  the  table  to  make  it  "complete. 


PART  IV:  0NED3P  COMPARISONS  WITH  AVAILABLE  SOLUTIONS 


21.  A  variety  of  ONED3P  calculations  were  performed  on  the  WES 
G-635  and  DPS/1  computer  systems  to  demonstrate  and  evaluate  how  well 
the  code  presently  functions.  These  demonstrations  are  presented 
herein. 

22.  Time  increments  and  loading  function  rise  times  which  led  to 
stable  calculations  were  selected  using  commonly  accepted  criteria;^ 
namely. 


At 
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max 
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max 
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where  Az  is  the  element  size  and  the  minimum  and  maximum  wave  speeds 
for  each  problem  are  determined  by  the  smallest  value  and  the  sum 

of  and  the  largest  value,  respectively.  Although  loading  func¬ 

tions  which  do  not  satisfy  Equation  12  may  be  used,  it  was  discovered 
that  such  calculations  generated  larger-than-expected  wave  speeds  for 
elastic  problems. 


Nonviscous  Problems 


Linear  elastic  single-layered  column 

23.  Figure  5  describes  the  problem  geometry,  material  properties, 
surface  loading  conditions,  and  boundary  conditions  for  four  0NED3P  test 
calculations  that  specify  linear  elastic  material  behavior.  The  first 
three  calculations  are  for  a  single  layer  of  material.  Problems  1  and  2 
were  computed  to  test  free  and  fixed  bottom  boundary  conditions,  respec¬ 
tively,  while  Problem  3  was  run  to  determine  whether  or  not  a  simple 
calculation  could  be  made  with  a  step  load  having  no  rise  time. 

24.  Stress-,  velocity-,  and  displacement-time  histories  for  each 
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Figure  5.  Problem  descriptions  for  linear  elastic  calculations 
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of  these  three  problems  are  contained  in  Figures  6  through  14.  Exact 
solutions  are  superimposed  on  selected  stress-time  histories  and  clearly 
show  that  0NED3P  can  handle  linear  elastic  calculations  very  well.  Elim¬ 
inating  the  loading  function  rise  time  in  Problem  3  resulted  in  slightly 
greater  stress  and  velocity  oscillations  than  those  in  Problem  2  as  well 
as  greater  wave  speeds  than  expected,  but  the  calculation  still  remained 
stable.  This  is  not  to  say  that  all  0NED3P  calculations  would  be  stable 
without  a  loading  function  rise  time.  Rather,  the  user  is  advised  to 

use  the  rise  timu  stability  criterion  (Equation  12)  for  all  computations. 

Linear  elastic  multilayered  column 

25.  Figure  5  also  contains  a  description  of  Problem  4,  which  is 

a  column  of  two  material  layers,  the  bottom  layer  having  twice  the  im¬ 
pedance  of  the  top  layer.  Results  for  this  calculation  are  shown  in 

Figures  15  through  17  and,  once  again,  the  code  does  an  excellent  job  of 

matching  the  analytical  solution. 

Linear  hysteretic 
single-layered  column 

26.  Salvadori  et  al.^  developed  an  analytical  solution  for  one¬ 
dimensional  stress  wave  propagation  through  linear  hysteretic  material 
which  was  applied  to  Problem  5  described  in  Figure  18.  Since  the  analyt¬ 
ical  solution  was  for  a  semi- inf inite  medium,  Problem  5  presented  an 
opportunity  to  test  the  transmitting  boundary  in  0NED3P  on  something 
other  than  a  linear  elastic  material. 

27.  Calculation  results  are  compared  with  analytical  results  in 
Figures  19  through  21.  Agreement  is  excellent.  Note  that  in  the  stress¬ 
time  histories  only  a  slight  bump  occurs  in  the  response  of  each  element 
at  the  times  when  waves  reflected  off  the  boundaries  would  normally  pass 
through  the  element. 

Nonlinear  hysteretic 
single- layered  column 

28.  The  only  available  solution  for  stress  wave  propagation 
through  highly  nonlinear  hysteretic  material  was  one  generated  by  the 
ONED  code  and  presented  in  Reference  8.  This  problem,  designated 


Figure  7.  Yalocity-time  histories  for  Problem  1 


qigp  rn  OISP  rM  OIGp 


Figure  8 


\n 


r 


dO!C  u  ;  JO  I G 


25 


DISP  -  CM  OISP  -  CM  DI9P 


0*0  25.0  50.0  U.O  100.0  125.0 

TIME  -  MSEC 


0.0  25.0  50.0  75.0  100.0  125.0 

TIME  -  MSEC 

Figure  14.  Displacement-time  histories  for  Problem  3 
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Figure  21.  Displacement-time  histories  for  Problem  5 
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Problem  b,  Is  described  in  Figure  22  which  also  shows  the  data  points 
used  in  both  codes  to  digitize  the  nonlinear  stress-strain  curve. 

29.  Stress-time  histories  from  the  two  code  calculations  are  com¬ 
pared  at  various  depths  in  Figure  23.  Obviously  the  first  compression 
wave  and  its  reflection  look  very  much  alike  using  either  code.  However 
it  appears  that  unloading  waves  which  result  in  low  stress  levels  travel 
much  faster  in  the  0NED3P  calculation  than  in  the  ONED  calculation.  The 
reason  for  this  lies  within  the  unloading-reloading  logic  used  by  each 
code . 

30.  Consider  Figure  24  which  shows  the  unloading  curves  for  any 
element  as  generated  within  each  code.  The  ONED  code  calculates  the 
stress  level  where  an  unloading  curve  bends  as  a  percentage  of  the  maxi¬ 
mum  previous  stress  computed  in  the  element.  Therefore,  if  the  unload¬ 
ing  curve  was  originally  defined  from  the  point  A,  then  an  unloading 
curve  from  point  B  would  look  like  that  shown  in  Figure  24a.  On  the 
other  hand,  the  0NED3P  was  designed  to  account  for  the  observation  that 
many  unloading  curves  bend  at  about  the  same  stress  level  regardless  of 
the  stress  value  from  which  they  originate.  The  resulting  unloading 
curve  from  point  B  as  computed  by  0NED3P  is  shown  in  Figure  24b.  Compar¬ 
ing  the  two  figures,  it  becomes  obvious  that  if  unloading  takes  place 
from  stress  level  B  to  stress  level  C  the  slope  of  the  unloading  curve 

at  C  in  the  0NED3P  calculation  would  be  greater  than  the  slope  at  C  in 
the  ONED  calculation.  Under  these  conditions,  unloading  waves  in  0NED3P 
would  travel  faster  than  similar  waves  in  ONED. 

Viscous  Problems 

Linear  viscoelastic  column 

9 

31.  Using  Laplace  transform  methods,  Morrison  solved  the  problem 
of  a  semi-infinite  column  of  linear  viscoelastic  material  subjected  to 

a  step  load  at  its  surface.  One  of  the  material  models  he  used  was  a 
three-parameter  model  like  that  of  Figure  la  where  the  three  parameters 
all  had  constant  values.  His  results  were  presented  in  a  nondimensional 
form  of  stress  versus  depth  in  the  column  at  constant  times.  By 


1.8288  M 


03  t  -  MSEC 


b)  SURFACE  LOADING  FUNCTIONS 


VERTICAL  STRAIN,  % 


c)  MATERIAL  PROPERTIES 

Figure  22.  Problem  6  description 
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cross-plotting  Morrison's  results  .it  given  depth.?  itlei  ini  erpol.it  ing 
additional  constant  Lime  curves  it  was  possible  l.u  generate  stress-time 
histories  I  rom  his  published  results  which  could  then  be  eomp.ircd  with 
0NKD3P  results  for  any  given  problem. 

32.  Figure  23  describes  the  linear  viscoelastic  0NED3P  calcula¬ 

tions  which  were  set  up  for  comparison  with  Morrison's  results  for  a 
given  set  ot  M ^  ,  M,  ,  and  ri  values.  The  analytical  results  assumed 
a  step  load  on  the  surface  of  the  column,  whereas  a  finite  rise  Lime  was 
chosen  for  the  0NED3P  calculations.  in  fact,  two  rise  Limes  were 

seLected:  1  msec  for  Problem  7  and  0.5  msec  for  Probiem  8. 

33.  Results  for  each  problem  are  shown  as  stress-time  histories 

at  various  depths  in  Figure  26.  Morrison's  solution  is  shown  as  dashed 
curves.  Although  the  finite  rise  times  in  each  problem  contribute  to 
poor  early-time  comparisons  at  each  depth,  the  late-time  comparisons 
look  very  good.  There  was  sufficient  viscosity  in  these  calculations 
to  cause  low-stress-level  waves  to  travel  at  a  speed  determined  by  the 
sum  of  and  and  to  arrive  at  each  depth  at  the  correct  time  (as 

predicted  by  Morrison's  solution) . 

3A.  Cutting  the  rise  time  in  Problem  8  cuased  an  overshoot  of 
calculated  stress  compared  to  Morrison's  predictions.  The  reason  for 
this  phenomenon  is  not  clear. 

35.  The  question  of  how  well  the  transmitting  boundary  works  in 
viscous  calculations  cannot  be  answered  in  this  section.  The  column  of 
material  in  Problems  7  and  8  is  long  enough  that  for  the  selected  mate¬ 
rial  a  stress  wave  would  not  reflect  back  from  the  bottom  boundary  to 
the  3-metre  depth  in  20  msec.  Transmitting  boundaries  for  viscous 
materials  will  be  discussed  later  in  paragraphs  44-46. 

Effect  of  loading  rate 
and  viscosity  on  wave  speeds 

36.  Consider  a  column  of  material  which  behaves  like  the  three- 
parameter  mechanical  model  shown  in  Figure  la  and  which  is  loaded  by  a 
loading  function  with  a  finite  rise  time.  Intuitively,  as  that  rise 
time  decreases,  stress  waves  travelling  through  the  column  should  travel 
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Figure  25.  Description  of  problems  for  comparison 
with  Morrison’s  solution 


faster.  The  reason  for  this  is  that  the  viscous  element  or  dashpot  be¬ 
lieves  more  and  more  like  a  rigid  element  as  the  rate  of  loading  on  it  in¬ 
creases  which,  in  turn,  magnifies  the  contribution  of  to  the  total 

stiffness  of  the  model.  As  this  happens,  the  slope  of  the  stress-strain 
curve  at  any  point  in  the  column  will  increase,  which  results  in  faster 
wave  speeds. 

17.  Now  consider  a  sinusoidal  loading  function  acting  on  a  three- 
parameter  material  having  constant  properties.  From  Kolsky^  the  speed 
ot  a  sinusoidal  dilations!  stress  wave  in  such  a  material  is  a  function 
of  its  frequent v  and  may  be  written 


(,  -  wave  speed 


2M1(M1  +  M2) 
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where  w  is  a  n>  .'mal  izto  frequency  and  is  equal  to  the  frequency  of 
tile  wave  ( -  —  t  )  times  5  he  ’  time  of  retardation"  of  a  Kelvin-Voigt  element 
(i)  which  in  this  case  is 


MiM2/<Ml  +  m2) 


Note  that  tor  is  i"  Low  I  re  queue  i  es , 
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waere  C  is  the  slowest  possible  wave  speed  in  the  material  and  is 
associated  with  tin  element's  quasi-static  behavior  whereas  for  very 

high  frequencies 
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where  C  is  an  upper  bound  on  the  wave  speed  and  is  associated  with 

ma  a  1 

the  parallel  spring  e  1 «  :i.ei,  l  s .  Finally,  combining  liquations  13,  14, 


and  I1),  one  has 
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F'te  Chat  for  a  given  set  of  and  M0  values,  C/C^  is  a  unique 

function  of  the  product  of  frequency  and  viscosity. 

38.  Equation  17  is  an  analytical  tool  which  can  be  used  to  pre¬ 
dict  tlie  speed  of  a  sinusoidal  stress  wave  in  a  three-parameter  material 
as  a  function  of  its  frequency  and  the  properties  of  the  material.  The 
question  is:  Will  stress  waves  calculated  by  0NED3P  behave  in  the  same 
way?  A  series  of  six  0NED3P  calculations  were  devised  to  answer  that 
question.  Those  calculations  are  described  in  Figure  27.  The  signifi¬ 
cance  of  the  parameter  values  and  loading  frequencies  which  were  chosen 
will  soon  be  apparent. 

39.  Two  examples  of  0NED3P  calculation  results  for  these  problems 
are  shown  in  Figure  28.  That  the  stress-time  histories  at  each  depth 

in  the  column  are  not  smooth  sinusoidal  curves  may  be  attributed  to 
numerical  approximations  inherent  in  any  finite  difference  or  finite 
element  code.  Observe,  also,  that  the  conditions  at  any  depth  do  not 
become  steady-state  until  after  approximately  two  stress  cycles.  All 
six  calculation  results  exhibited  similar  behavior. 

40.  Focusing  attention  on  the  stress-time  histories  at  the  three 

greatest  depths,  the  following  method  was  applied  to  determine  the  wave 
speed,  t  ,  for  each  problem.  Utilizing  the  common  drafting  technique 
for  drawing  parallel  lines  with  two  triangles,  a  line  was  chosen  for 
each  problem  which  best  described  the  intersection  of  each  stress  cycle 
with  the  time  axes  at  the  three  greatest  depths.  Naturally  this  was 
done  only  for  stress  cycles  which  occurred  after  the  stress  wave  became 
steady.  The  inverse  slope  of  this  line,  being  the  speed  with  which  each 
stress  cycle  propagates  along  the  column,  was  then  divided  by  from 

Equation  15  and  the  results  were  plotted  in  Figure  29  which  contains  the 
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Figure  27.  Description  of  problems  used  to  study  the  effects  of  fre¬ 
quency  of  a  sinusoidal  stress  wave  on  its  wave  speed 
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unique  curve  determined  by  substituting  the  and  M2  values  for 

these  six  problems  into  Equation  17. 

41.  Qualitatively,  the  results  of  these  six  calculations  were 
quite  good  in  that  they  form  points  on  a  curve  which,  considering  the 
nature  of  the  0NED3P  solution  algorithm,  ;Ls  unique  and  which  is  nearly 
parallel  with  the  predicted  curve.  Quantitatively,  however,  the  calcu¬ 
lated  response  is  shifted  to  the  right  of  the  predicted  response.  This 
may  be  looked  at  in  one  of  two  ways:  first,  for  a  given  frequency  and 
viscosity,  the  calculated  speed  of  a  sinusoidal  stress  wave  is  less  than 
would  be  predicted  by  theory,  or,  second,  in  order  for  a  sinusoidal 
stress  wave  to  travel  at  a  given  speed,  it  either  must  have  a  greater 
frequency  than  theory  would  require  or  it  must  be  in  a  material  which 

is  more  viscous  than  theory  would  dictate  (or  both).  Time  and  funding 
do  not  permit  a  more  thorough  examination  of  these  results.  However,  it 
is  suggested  that  if  the  discrepancy  between  computed  and  predicted  wave 
speeds  is  due  to  numerical  approximations  within  the  code,  a  calculation 
with  a  much  finer  grid  and  time  step  (and  therefore  more  expensive)  might 
result  in  a  better  correlation. 

42.  Although  sinusoidal  loading  functions  are  analytically  clean 
in  the  sense  that  they  possess  only  one  frequency  component,  the  types 
of  dynamic  loading  functions  which  are  used  in  laboratory  testing  or 
which  are  observed  in  field  experiments  contain  an  infinite  number  of 
frequency  components  and  may  be  approximated  in  many  cases  by  a  step 
load  with  a  finite  rise  time  such  as  that  shown  in  Figure  30b.  As  a 
further  exercise  in  studying  the  effects  of  viscosity  on  wave  speed,  the 
three  problems  described  in  Figure  30  were  calculated  and  the  resulting 
wave  forms  were  plotted  in  Figure  31.  Clearly,  for  a  given  rise  time, 
increased  viscosity  results  in  faster  wave  speeds.  The  reason  for  this 
is  shown  in  Figure  32  where  it  can  be  seen  that  increased  viscosity 
causes  a  stiffer  stress-strain  response  in  a  given  element.  If  one  took, 
the  rise  time  of  the  loading  function  as  one-fourth  of  the  period  of  a 
sine  wave,  it  would  be  possible  to  calculate  C/Cg  values  to  be  plotted 
in  Figure  29.  However,  the  initial  arrival  times  for  sinusoidal  stress 
waves  did  not  correlate  well  with  theoretical  predictions;  neither  do 
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Figure  30.  Description  of  problems  used  to  study  the  effects  of 

viscosity  on  wave  speed 
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the  arrival  times  for  Problems  15,  16,  and  17. 

4S.  At  this  time  it  is  not  recommended  that  Equation  17  be  used 
to  predict  initial  arrival  times  for  stress  waves  generated  by  loading 
functions  like  that  in  Figure  30b.  It  is  felt  that  the  infinite  fre¬ 
quency  components  in  a  finite  rise  time  step  load  violate  the  assumptions 
under  which  Equation  17  was  derived.  However,  Figure  31  does  show 
clearly  that  wave  speeds  in  a  dynamic  problem  with  a  finite  rise 
time  are  a  function  of  viscosity.  Furthermore,  it  will  be  demonstrated 
in  Part  V  that,  for  a  given  viscosity,  loading  rates  on  the  order  of  a 
fraction  of  a  millisecond  can  affect  wave  speeds. 

More  on  the  transmitting  boundary 

44.  The  0NED3P  transmitting  boundary  was  shown  to  work  quite  well 
for  the  nonviscous  linear  hysteretic  problem — Problem  5.  As  a  check  on 
how  well  it  handles  viscous  materials.  Problem  5  was  recalculated  with 
an  M.;  value  equal  to  69  MPa  and  three  different  n  values  (see  Fig¬ 
ure  18  for  other  problem  parameters).  These  calculations  were  assigned 
n  values  equal  to  10,  100,  and  1000  MPa-msec,  respectively. 

45.  Results,  in  the  form  of  stress-time  histories,  are  shown  in 
Figures  33,  34,  and  35.  Assuming  that  the  stress-time  history  at  any 
depth  should  be  a  smoothly  decaying  function  after  the  peak  stress  has 
been  reached,  it  is  obvious  that  some  energy  is  reflected  from  the  bot¬ 
tom  boindary.  For  the  highest  viscosity  the  transmitting  boundary 
appears  to  behave  too  stiff,  resulting  in  a  small  compression  wave  being 
reflected  back  to  the  top  of  the  column. 

46.  It  is  left  to  the  discretion  of  the  0NED3P  user  whether  or 
not  to  use  the  transmitting  boundary  for  his  problem.  These  three 
calculations  are  only  offered  as  examples,  and  it  would  be  improper  to 
draw  from  them  general  conclusions  concerning  all  viscous  problems. 
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Figure  33.  Stress-time  histories  for  Problem  5; 
n  =  10  MPa-rasec,  M„  =  10,000  MPa 


PART  V:  0NED3P  SIMULATION  OF  FIELD  EXPERIMENTS 


Background 

47.  Several  field  experiments  have  been  conducted  by  WES  person¬ 
nel  to  test  the  response  of  shallow-buried,  flat-roofed,  concrete  box 
structures  to  an  approximate  plane  wave  surface  load  generated  by  soil 
berm-covered  explosions.  The  geometries  of  these  experiments  are  pre¬ 
sented  in  References  2  and  3.  Soil  stress  histories  at  different  depths 
were  measured  in  the  field  during  the  events. 

48.  Laboratory  material  property  tests  were  conducted  on  the 
sand  backfill  placed  around  and  above  these  structures  to  determine  the 
sand's  compressibility.  Rapid  loading  of  sand  backfill  laboratory 
samples  resulted  in  stiffer  stress-strain  responses  than  those  obtained 
with  static  loadings.  Field  test  results  showed  that  the  downward 
propagating  stress  waves  travelled  through  the  backfill  sand  with  a 
speed  much  faster  than  the  p-wave  velocities  determined  from  seismic 
refraction  surveys.  It  has  been  suggested  that  this  phenomenon  can 
also  be  explained  by  the  rate-dependent  behavior  of  the  sand. 

49.  An  effort  was  previously  made  to  simulate  these  field  experi¬ 
ments  with  the  ONED  code  using  different  rate- independent  stress-strain 
curves  to  simulate  the  rate-dependent  behavior  of  the  sand  backfill  at 
different  depths. ^  This  brute-force  method  of  accounting  for  rate 
dependence  is  at  best  an  art  which  requires  some  posttest  knowledge  and 
is  therefore  not  useful  for  pretest  predictions.  It  would  be  both  more 
convenient  and  physically  more  sound  to  use  a  true  rate-dependent  one¬ 
dimensional  wave  propagation  code  to  try  to  simulate  these  field  experi¬ 
ments.  In  this  regard  the  development  of  0NED3P  is  very  timely. 

50.  To  demonstrate  how  0NED3P  can  be  applied  to  real  problems, 
backfill  response  in  two  of  the  above-mentioned  field  experiments,  desig¬ 
nated  as  FH4  and  FH5,  was  simulated  with  the  new  code.  Backfill  proper¬ 
ties  were  assumed  to  be  the  same  in  both  experiments  and  only  the  sur¬ 
face  loading  functions  were  different. 
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Evaluation  of  Model  Parameters 


51.  At  present,  the  most  sound  method  for  determining  quantita¬ 
tive  values  for  the  three  mechanical  parameters  required  by  the  0NED3P 
code  is  to  select  those  parameters  which  will  best  reproduce  both 
"static"  and  "dynamic"  laboratory  uniaxial  strain  test  data.  This  can 
be  done  by  using  a  simplified  version  of  0NED3P  called  a  driver.  Put 
simply,  the  driver  determines  the  strain  response  of  the  three-element 
model  due  to  a  known  applied  stress-time  history. 

52.  Consider  the  0NED3P  mechanical  model  and  typical  labora¬ 
tory  data  shown  in  Figure  36.  In  order  to  properly  evaluate  the  model 
parameters,  one  must  have  a  "quasi-static"  response  curve  and  at 

least  two  (but  preferably  more)  "dynamic"  stress-strain  curves  that  will 
clearly  depict  the  effects  of  rate  of  loading  on  the  stress-strain  re¬ 
sponse.  The  "dynamic"  test  data  should  include  the  fastest  loading 
rates  expected  under  field  conditions.  The  "quasi-static"  data  can  be 
derived  from  much  slower  loading  rates;  e.g.,  rise  times  on  the  order 
of  tens  of  seconds  are  usually  adequate.  The  three  mechanical  parameters 
will  now  be  examined  one  at  a  time. 

53.  First  there  is  the  function.  It  is  defined  by  the 

"quasi-static"  curve  since  very  slow  loading  rates  effectively  cause 
the  viscous  and  M2  spring  elements  in  Figure  36  to  disappear. 

Since  0NED3P  is  an  incremental  code,  M^  is,  in  fact,  the  tangent  slope 
of  the  quasi-static  curve  at  any  point.  Thus,  for  some  nonlinear 
hysteretic  materials,  may  vary  dramatically. 

54.  Next  is  the  value  of  M2  ,  which  in  the  current  version  of 
0NED3P  is  taken  to  be  a  constant.  For  extremely  rapid  loading  rates  the 
viscous  element  is  nearly  rigid  and  the  response  of  the  three-parameter 
model  is  governed  by  the  sum  of  M^  and  M2  .  Thus,  one  possible  means 
of  establishing  M2  is  to  determine  the  value  of  the  initial  tangent 
slope  of  the  stiffest  dynamic  stress-strain  loading  curve  and  call  that 
Mf  +  M2  .  This  tangent  should  form  an  upper  bound  to  all  of  the  avail¬ 
able  dynamic  stress-strain  data.  M2  can  then  be  calculated  because  the 
value  of  M^  at  the  origin  is  known. 
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n  ,  which, 


55.  Finally  there  is  the  viscosity  coefficient, 
again,  is  presently  taken  to  be  a  constant  in  0NED3P.  Since  n  is 
effectively  zero  under  quasi-static  loading  conditions  and  is  effec¬ 
tively  equal  to  «>  under  almost  instantaneous  loading  conditions, 

the  task  is  to  find  the  most  appropriate  value  of  n  (between  zero  and 
")  that  will  produce  agreement  with  the  available  "dynamic"  stress- 
strain  data.  This  can  be  done  by  a  tr ial-and-error  process  (once  M^ 
and  are  specified).  That  process  entails  varying  n  and  executing 

the  driver  with  each  dynamic  forcing  function  associated  with  the  labora¬ 
tory  tests. 

56.  For  example,  consider  the  laboratory  uniaxial  strain  test 
data  derived  from  tests  on  samples  of  sand  backfill  used  in  the  FH4 
field  experiment.  These  data  are  shown  in  Figure  37.  The  "static"  re¬ 
sponse  of  the  FH4  sand  (Figure  37a)  is  well  behaved,  and  a  bilinear  ap¬ 
proximation  to  the  average  of  the  loading  curves  was  assumed  for 
while  the  unloading  response  was  assumed  to  be  a  straight  line.  The 
change  in  slope  of  the  loading  curve  occurred  at  20  MPa  and  7.6  percent 
strain,  while  the  unloading  curve  was  drawn  from  35.2  MPa  and  11.8  per¬ 
cent  strain  to  0  MPa  and  10.8  percent  strain. 

57.  As  for  M;  ,  the  steep  tangent  shown  in  Figure  37b  was  chosen 
at  the  upper  bound  to  the  data;  its  slope  was  such  that  M^  was  calcu¬ 
lated  to  be  approximately  10,000  MPa. 

58.  Finally,  several  values  of  n  were  arbitrarily  chosen  and 
the  driver  was  exercised  with  the  applied  load  functions  from  laboratory 
tests  D4.6,  D4 . 7 ,  and  04. 8,  which  are  all  shown  in  Figure  37b.  A  value 
of  n  =  10  MPa-msec  was  finally  selected.  It  gave  a  good  approximation 
to  the  D4.7  dynamic  stress-strain  curve  but  significantly  undercut  the 
initial  slopes  of  the  faster  tests  (D4.6  and  D4.8)  as  shown  in  Figure  38. 
Before  discussing  the  0NED3P  results  using  these  parameter  values  it 
should  be  noted  that  another  possible  method  for  establishing  the  mate¬ 
rial  parameters  is  as  follows.  Given  some  M^  or  "quasi-static"  func¬ 
tion,  M^  and  n  could  be  derived  through  iterative  0NED3P  calculations 
against  experiments  like  FH4  such  that  the  arrival  times  of  the  stress 
wave  fronts  eventually  match  those  measured  in  the  field  tests.  Other 
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Figure  37.  Laboratory  uniaxial  strain  test  results  for  FH4  backfill  sand 
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Figure  38.  A  comparison  of  the  three-parameter  model 
fit  with  dynamic  FH4  laboratory  results 
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specific  problem  needs  not  covered  herein  could  dictate  additional  ap¬ 
proaches  to  determining  the  0NED3P  parameters. 


Calculation  Results 


59.  Two  0NED3P  calculations  were  set  up  for  comparison  with  FH4 

3 

and  FH5  field  test  results.  Sand  density  was  taken  as  1685.3  kg/m  . 

A  grid  spacing  of  0.3175  cm  was  used  and  the  time  increment  was  taken 
to  be  0.001  msec.  In  fact,  the  problem  set-ups  were  identical  to  the 
previously  mentioned  ONED  code  simulations  of  the  same  tests.  (0NED3P 
ran  30  percent  faster  than  ONED  for  each  calculation.) 

60.  Figures  39  and  40  compare  stress-time  histories  generated  at 
various  depths  by  0NED3P  with  the  measured  field  stress  wave  forms. 

Also  included  on  those  figures  are  the  dynamic  stress-strain  curves 
generated  by  0NED3P  at  the  corresponding  depths.  Although  these  simula¬ 
tions  are  only  reported  as  an  example  of  how  to  apply  0NED3P  to  a  real 
problem,  a  brief  discussion  of  these  results  is  still  warranted. 

61.  First  it  is  obvious  that  the  calculated  wave  forms  are  differ¬ 
ent  than  the  measured  wave  forms.  There  are  at  least  two  things  which 
could  be  done  to  improve  that  comparison.  One  would  be  to  obtain  a  bet¬ 
ter  fit  to  the  dynamic  laboratory  data  in  the  previous  section  using  the 
model  in  its  present  form.  However,  there  is  presently  not  enough  flex¬ 
ibility  in  the  model  parameters  to  preserve  both  the  low-  and  high- 
stress  stress-strain  responses  shown  in  Figure  38.  On  the  other  hand,  it 
is  very  likely  that  further  developments  in  0NED3P  might  result  in  a 
better  simulation  of  the  measured  wave  forms.  For  example,  if  viscosity 
increased  with  stress  level,  a  sharper  wave  front  would  be  calculated. 
Such  modifications  to  0NED3P  are  being  considered  but  are  not  within  the 
scope  of  this  report. 

62.  On  a  more  positive  note,  the  stress  wave  arrival  times  at 
various  depths  were  calculated  with  greater  success.  Note  that  the  cal¬ 
culated  stress  waves  slow  down  as  they  travel  deeper  into  soil.  The 
observation  is  consistent  with  the  dynamic  stress-strain  behavior  of  each 
element.  As  the  wave  travels  deeper  into  the  soil,  the  dynamic 
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stress-strain  response  approaches  the  "static"  response.  Initial  stress 
wave  velocities  predicted  by  quasi-static  laboratory  data  would  have 
been  much  smaller  than  those  measured  in  the  field. 

6J.  Note  also  that  the  wave  speeds  calculated  for  the  FH5  field 
test  simulation  are  greater  than  those  for  the  FH4  simulations. 

Since  the  only  difference  between  the  two  calculations  is  the  rate  of 
loading  (40  microseconds  to  21.3  MPa  for  FH4  and  30  microseconds  to 
127.7  MPa  for  FH5) ,  the  results  demonstrate  conclusively  that  the  rate  of 
loading  must  influence  wave  speeds  by  causing  a  stiffer  dynamic  stress- 
strain  in  the  material. 


PART  VI:  SUMMARY 


64.  Recent  experimental  data  on  soils  showing  that  the  application 
of  loads  with  submillisecond  rise  times  results  in  significant  rate- 
dependent  compressibility  behavior  have  prompted  the  need  for  a  one- 
dimensional  plane  wave  propagation  code  for  layered  rate-dependent 
hysteretic  or  visco-compacting  materials.  Accordingly  an  explicit  one¬ 
dimensional  finite  element  code  named  0NED3P  has  been  developed  which 
incorporates  a  three-parameter  (spring  and  dashpot)  mechanical  model 
consisting  of  a  linear  spring  and  dashpot  in  series  coupled  in  parallel 
with  a  piecewise  linear  hysteretic  or  compacting  spring. 

65.  0NED3P  allows  for  the  analysis  of  a  column  of  multilayered 
visco-compacting  soils  loaded  by  a  digitized  surface  pressure-time 
history.  Any  set  of  consistent  units  may  be  used  with  this  code  and  re¬ 
sults  may  be  obtained  in  the  form  of  stress-,  strain-,  acceleration-, 
velocity-,  and/or  displacement-time  histories  as  well  as  stress-strain 
curves  using  standard  Calcomp  software. 

66.  Several  demonstration  problems  were  calculated  using  0NED3P 
to  evaluate  its  features  and  capabilities  against  (a)  available  analyti¬ 
cal  solutions  for  viscous  and  nonviscous  problems,  (b)  other  code  solu¬ 
tions,  and  (c)  measurements  from  field  experiments  that  evince  rate- 
dependent  soil  behavior.  In  general,  results  were  extremely  good. 

67.  Special  attention  was  given  to  the  effects  of  loading  rate  or 
frequency  on  wave  speeds  in  viscous  media  and  to  methods  of  deriving 
the  0NED3P  model  parameters  from  laboratory  material  property  test  data. 
Program  listings  and  a  user's  guide  for  0NED3P  are  included  in  Appendix  A. 
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APPENDIX  A:  WES  USER  GUIDELINES  AND  PROGRAM/DR 1VER  LISTINGS 


ON ED 3 P 


1 .  ONED3P  is  set  up  for  running  on  the  WES  GE  635  computer  system 
CARDIN  subsystem*  and  requires  a  free-field  data  file  in  the  form  of  a 
permanent  file  for  execution,  whose  name  should  appear  on  line  number 
8030  of  0NED3P : 


8020$ :PRMFL: 39, Q,L,ROSDJOC/ file  name 

If  there  is  more  than  one  input  quantity  on  a  line  of  the  data  file 
those  quantities  should  he  separated  by  commas  or  blanks.  Any  line  num¬ 
bers  used  to  generate  the  file  must  be  stripped  before  0NED3P  can  be 
executed.  Hollerith-type  information  must  be  enclosed  in  double  quotes. 
The  format  of  a  typical  data  file  follows. 

Section  1:  General  Information 

Line  "Problem  Title" 

Line  Boundary  condition 

The  boundary  condition  is  an  integer  which  selects  the  type  of  bottom 
boundary  condition  and  may  have  the  following  values: 

1  =  free  2  =  fixed  3  =  transmitting 

Section  2:  Material  Properties 

"Title  Describing  Materials" 

Number  of  layers  in  the  soil  column 

*  Execution  time  for  any  calculation  on  the  GE  635  may  be  figured  by 
assuming  that  ONED3P  needs  0.00007  hundredths  of  an  hour  of  CP  time/ 
element/time  step. 


Line 

Line 


A1 


Line 

Layer  number,  number  of  elements  in  layer, 
density,  layer  height 

,  layer 

Line 

n  ,  M0  for  that  layer 

Line 

Number  of  stress  and  strain  pairs  (of  o-G 

coordinates) 

defining  the  piecewise  linear  segments  of  the  Mj 
load  curve  including  the  origin  number  of  stress  and 
strain  pairs  defining  the  unloading  curve 


Line(s)  Values  of  stress-strain  pairs  defining  load  curve 

in  sequential  order  beginning  at  the  origin 

Line(s)  Values  of  stress-strain  pairs  defining  Mj  unloading 
curve  beginning  at  user-chosen  point  on  the  load  curve 


The  last  five  (or  more)  lines  are  repeated  until  all  of  the  layers  are 
defined  beginning  with  layer  1  at  the  top  of  the  soil  column. 


Section  3:  Surface  Loading  Function 

Line  "Title" 

Line  Number  of  stress-time  data  pairs 

Line(s)  Stress-time  data  pairs  beginning  at  the  origin 

Section  4:  Execution  and  Print/Plot  Parameters 

Line  Problem  time  at  which  calculation  will  stop,  time 

increment 

Line  Print  interval,*  plot  interval,*  number  of  locations  in 

the  column  for  which  plot  data  are  to  be  saved** 

Line  Numbers  of  the  elements  or  nodes  in  the  column  for 

which  plot  data  are  being  savedt 

Line  Total  number  of  plots  requested 

Line  Type  of  plot,  element  or  node  number 


*  Print/plot  interval  takes  an  integer;  1  means  save  data  at  every  time 
step  for  print/plot;  4  means  save  data  every  fourth  time  step,  etc. 

**  The  dimension  of  the  FPLOT  vector  in  0NED3P  (line  175)  limits  the 
amount  of  plot  data  which  can  be  saved.  That  dimension  must  be 
greater  than  500  plus  six  times  the  number  of  time  points  to  be  saved 
times  the  number  of  elements  being  saved. 

+  For  each  number  listed  results  are  saved  for  both  the  element  and 
node  having  that  number. 
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Line  X  value  of  plot  origin  in  inches,  Y  value  of  plot 

origin  in  inches,  plot  angle  in  degrees,  X  axis  scale 
factor  in  units/inch  of  plot,  Y  axis  scale  factor, 
length  of  X  axis  in  inches,  length  of  Y  axis 


The  last  two  lines  of  data  are  repeated  sequentially  until  all  requested 
plots  have  been  described.  These  lines  require  the  following  explana¬ 
tion.  First  of  all  the  plot  type  is  an  integer  with  the  following 
values : 


Value 

1 

2 

3 

4 

5 

6 


_ Type  of  Plot _ 

Stress-time  history 
Strain-time  history 
Stress-strain  curve 
Acceleration-time  history 
Velocity-time  history 
Displacement-time  history 


Next  the  origin  of  each  plot  refers  to  its  position  on  a  standard  34-in. 
Calcomp  drum  plot  where  the  origin  on  the  drum  plot  is  located  near  the 
bottom  edge  of  the  paper.  The  plot  angle  may  be  either  zero  degrees  for 
plots  whose  X  axis  is  parallel  to  the  bottom  edge  of  the  drum  roll  or  90, 
180,  or  270  degrees  rotated  in  a  counterclockwise  direction. 

2.  As  an  example  of  a  typical  data  file  for  0NED3P,  consider 
Figure  A1  which  shows  the  file  that  was  used  for  one  problem  in  Part  IV. 

3.  A  listing  of  0NED3P  follows. 
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"MOP F  I  I- Oh  CHECK  MO.  1  ■'  j 
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c!  5  c! 
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3  I 
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1  i-  5  0 

d .  y  15.  y  U .  y  4 .  y  1  1.1,  y  5  •  »  1  . 

<S 

Figure  Al.  ONED3P  input  file  for  problem  7, 
Part  IV 


SECTION  1 


SECTION  2 


SECTION  3 


SECTION  4 


*  -  i.  ' :  % 


»*«***■  «r4W<fc,  v«4  ./ 


due-  t'  :p 


file  fhi.iE  no 


1 


1  II :  5:  *:  l  IJF:  66  6  1 


"01 n  1  1  :  <  T 

imli'f:  I  DENT :  f  ii .  I>  JQC  <  CNP'T  I  L 
""1  4  1  :  U !  EF  I  D:  P"  D  JQf  lx  MPT  I : 
i.m  It- i:  OFT  ION:  FOPTPhM 
""1  31 :  F OF'T'i' : :  FEF 
in 1 1  "1. :  L  1 M  I  T  ; :  1  n.  5>.l 

i.i  1 i^'  in. 

On  t.*C  FPOC’FhM  ONE  DOF' 

1 1 1 14  III. 

H i if. m;  h  ome-dimem: iomhl  vi  :i;oelm:th:  mhve  pfopfntht i on  lode 

iiiiimu;  THL  (DH'llTMTlVE  BEHhVIOF  mhY  BE  FEF'FE I EMTEB  BY  ft  MECHANICAL 

0  n  7  Oi.  MO  BEL  COM-I.TINC*  OF  H  FhPhLLEL  IFF INS  HMD  MAXWELL  ELEMENT 

ii  0;;.  iii' 

ii 4 in:  limit:  hpe  amv  e t  of  con: i stent  limit::  fof  instance* 

*.n »•»>*;  mehhf-h . r ML :  •  m:ec •  metef: .  m3 

"11  ML 

ii  1  ill  i  HhF  hi.  TEP*  1  AN  :•  TITLED*  FI  LE*64<  F  T  I  TL*40 

" i  :  ii  r armor <  :l  •  1  o* 50 • « : u •  1  o.  so ■  •  el  •  1  o»50 ■  > eu •  1 0*50 •  •  eme  < 65  o<  1 0 > 

II 1  31  .  IFLhh  .  65  0  .  .  LM •  M  .L  •  1  0  •  .  m:U  •  1  0  • 

014  0  BIMEM:  I  OH  TIME  » 700  •  •  BUN!  •  700-* .  BUM6  <  700  •  *  I  S  *  4  •  .EPS  <  3«  650*  • 

11 1 45  I  IS' 3.650'. 

"150  ftC  1;  *E<650.  <  VEL  •  c’<  *50  *  <I'I  iP'6.650* .AMASS <650'  «HLftY<650*  <SMAX<65U< 

i.i  16  0  K  tMH  ’  *  650  '  <  Ml  EEP  *  65"  .<  <LS  1 611  • .650  1  <  Z  <  650  •  »ZS  <  650 »  <ETft<  1  O  '  <  HM2  <  10*  < 

"1  ro  :■  :  TFE  ■  3  I  <  :  TPh  ■  3  *  <F  <  65 0  <  <  SMI  <  650  '<1110' 650  ' 

0 1  75  B 1  MEM :  I  DM  FF'LOT  < 65  00 0  < 

01*0  CALL  F:  :OF'T  I  67<  1  <  1  <  0  • 

"6  i.i  OC 

0,3  05  IFL'3B=0 

i  1 6  F  "  FEhB  loEO1  PTITL 

"6*0  F'FIMT. 

ii  3 ii 0  F'FIMT* 

n  1  n  F'FIMT. PTITL 

ii  36  0  F'FIMT. 

1 1  3  3 OL  F'EhB  BOTTOM  BOUNDARY  CDI'E 

0336L  BOUNDARY:  1=FFEE<  6=FIXED.  3=TPftNSMITTINS 

0  34  0  F'EftB  ■  3*<  1  06  0  ■  MB  TYPE 

0  346C 

0  34  3C  :ET  :TF'hIH  INCREMENT  TOLERANCE 

0  3440 

0  345  TDEF'7  =  6.E-6 

0 346  F'F'IMT.  "  7 T Fft I N  INCREMENT  TOLERANCE*"  <  TBEPS 

1 1 350  BO  10  I  =  l<3 

0 3*  i.i  BO  10  .1=3  <  F. 5  ii 

0  3  7  o  :  1 6  *  I .  J  '  =  0 . 

0330  EF;  *  I  .  J  •  =  0. 

0  3*0  10  FONT  I HUH 


u4  0  M 

BO  11  1=1.6 

n4  t  i> 

BO  11  J=1.6 

» 1 4 1 1 

ftC  C  •  I .  J  *  *  0 . 

i  »4  7:  • 

VEL  ■  I .  J* =0. 

.14411 

b  i :  p  •  i .  ,i  • = o 

1(4*7- «  ' 

11  CONTINUE 

ll4r.  H 

BO  16  1=1.6 

m4  7h 

F*  I  ■=  0. 

i(4  U 

ML  H  i  ■  I  i  =11 

A5 


□Hfc  li -’P 


COMT 


FILE  PHbB  HD 


048  i.i 
ii*.  ii  i.i 
i  *5  l  0 

n  sj:  o 
ft  5  5  i’i 
054  0 
0544 
U54* 
055  or 
OS6  0»I 
ill:. }:  Of 

064  nC 
0r-5»  UC 

0*6 >.i 

ft*.  "*  0 

i.lr.X  I.I 
0*80 
i.i  TOO 

on  oc 

075  OC 
nr?ui; 
074  OC 
0750 
0760 
0  7  SO 
0780 
Ox  0  0 
Ox  1  n 
Ox  5  0 
08  30 
0x4  0 
0850 
0x6  0 
0870 
i.i  xx  0 
1 1 8  8  0 
i.i  8  0  0 
081  0 
0850 
08  ?:  0 
084  0 
0850 
0860 
087  0 
088  0 
0880 
1  000 
1  0 1  0 
l  o  y.  o 

1  04  o 
1 05  o 
l  o*  o 
1  07  0 
1  Oxo 
t  080 


10:  s: :i  08  86  81 


hMh ' T • 1 • =  0 . 
:mh:  •  i  •  =o. 

: mi • i ■ =  o . 
15  EMh: -:»I»=n. 

E  1  =  0. 

E5=  0. 

HP1  =  ft. 

HP 5*  0. 


COMPUTE  LUMPED  MASS  VALUES  AMD  INPUT  LAYER  MATERIAL  PPOPEPTIE 

READ  1  S?>  l  0S0>  TITLE 
PRINT. TITLE 
READ  '  S9«  1  OSO'  NL 
nma ; : = i 
DTMIN* 1 .  E5 


1_M=  LAYER  NO.!  NZ=  NO.  DF  ELEMENTS  IN  LAYER 
P0=  LAYER  DENSITY!  H=  LAYER  THICKNESS 

READ  ■  S  r*.  1  Oc 0 >  LN .  NZ .  RO >  H 

IF'LN.NE.l'  STOP 

READ  '  ;•?.  1  OSO  ■  ETA  •  LN  ■  .  AMS  -  LN> 

READ'S'?.  1  OS  0  >  NSL  '  1  '  .NS  U  -  1  > 

READ  •  39 .  1  OS  0  •  •  S L  •  1 .  I  >  « EL  •  1 .  I  '  .  I  =  1  *  N S L  • 1  >  > 

READ  ■  39 .  1  Oc'  O'  •  S  U  •  1 .  I  '  .  EU  1  •  I  >  .  I  =  1  •  NS  U  •  1  >  > 

DZ=  H  FLOAT 'NZ' 

ama:=  pq*dz  s. 

: LOPE*  '  S U • 1 . 1 ■ - SU ■ 1 . S > > * ■ EU ■ 1 . 1 > -EU '. !»£>>♦  AMS < 1 > 

CP=  SORT- SLOPE  RO' 

DT  =  DZ  CP 

IF ' DT. LT . DTMIN 1  DTMIN*DT 
Z  ■  1  '  =  0 .  0 

IF  •  NZ .  E 0 .  0  ■  i>0  TO  1* 

DO  14  1=1. NZ 

AMA S : ' I ' =  AMh 7  S  ' I  ' +  ANA  S 

ama : : ■ i ♦ i ■ =  ama : s • i + i ' +  amas 
Z  ■  I  + 1 ' =  Z' I ' +  DZ 

z: '  i  ■=  z- 1  '+  dz  s. 

14  ML  A  Y ' I • =  LN 
nma : : =  i + i 
NEL=  NZ 

1C.  IF  NL.EC'.l'  >S0  TO  cl 
DO  c'O  J=S.HL 

READ • S?. 1  OS  0 ■  LN.NZ.PO.H 

IF'LN.NF.J'  STOP  "LAYER  INFO  OUT  OF  ORDER" 

READ • S?« 1  OS 0 •  ETA ■ LN • . AMS ' LN • 

READ • S?«  1  OS 0 •  N S  L • LN  > • NS U ' LN > 

READ ' S?. 1  OS 0-  • SL ' LN. K • «  EL ' LN-K • . K  =  1 > NSL ’ LN • > 

READ ■ S?> 1  OS O'  • SU'LN.K ' » EU • LN. K  ■ »K  =  1  -NSU 'LN' > 

DZ=  H  FLOAT 'NZ' 

AMA I =  RO*IC  S. 

:lofe=  ■  :u  •  ln.  i  •  -:u  ■  ln.s  •  •  ■  eu  •  ln.  i  •  -eu •  ln«s>  »  +ams  <ln.< 


A6 


□HE  IMP 


I'.OHT 


1  (IS  5:31  (IS  £'3  y  j 


FILE  PAGE  HD. 


1  1  0 1 1 

CP=  I  OPT  '  ;  LDF'E  P'0  ’ 

1  1  1  0 

dt=  rccp 

1 1  £  0 

IF'DT.LT.DTMIN'  DTMI  N=DT 

1  1  30 

no  13  i=i. re 

11411 

ana:  s  iNMfi : :  •  =  ama  : :  •  hma  : s  >  +  ama- 

1  I5u 

AMAS  S  >  NMAS  S  +  1  •  =  HHA  "  3  •  NMAS  3  +  1  ■  +  AHA  ' 

1  1  *  u 

£  •  r imh  : + 1 .  =  z  •  hmh  s :  ■  ♦ic 

1  1  7  0 

z :  •  nmass  >  =  z  ■  hmh:  s  >  +  nz*o.  5 

1  13  ii 

hlay-nma: : •  =  lm 

1  iso 

nma::  =  nma  r 7+1 

1  c!  U  M 

13 

COHT 1 HUE 

1  £  1  0 

HEL=  HEL+  HZ 

1  ££  0 

c‘0 

COHT I HUE 

i  £ o 

£1 

IF  ‘MISTYPE.  EC'.  3 >  AMASS  ■  NMASS. '  =  AMASS  .  NMASS  >  +  AMR' 

1  £4  0 

C  P I H I T=C P 

1  £5  0 

L  F'LA :  T=  CP 

1  M  0 

PPIHT  1  0  30.  '  I  •  Z  '  I  '  «  AMASS  ■  I  >  .  1  =  1 .  NMA" :• 

1 

IWJM1  <  1  •  =0. 

1  3  30 

DUM£  •  1  ■  =  0 . 

1  34  0 

NF'LT=£ 

1  .:5  ML 

1  33  OC 

INPUT  SURFACE  STRESS  TIME  HISTORY 

i  s  roc 

1  330 

THYME=  0. 0 

1  330 

F ' 1 1 =  S  STPS  <THYME. 1> 

140  DC 

1 4 1  OC 

ES  TAEL  I SH  A  TABLE  OF  MOIUJL I  FOP  THE  Ml  FUNCTION 

1 4c  OC 

1 4  3  0 

I'D  £3  1  =  1.  ML 

1 45  0 

no  £4  N=£ ,  NSL  • I » 

1430 

; LP=  ■  S  L  *  I « N  *  - SL • I . M- 1 ■ »  •  'EL ■ 1 , N  >  -EL  Cl ,  N- 1  » , 

147  0 

IF  <SLP.LT . 0. ’  STOP  "  NEGATIVE  LOADING  MODULUS" 

1 43n 

£4 

:l'I,h-i»=  :lp 

1430 

HM1=  NSL'I’-l 

15  30 

DO  £3  N=£ . HSU  1 I 1 

1 540 

SLP=  <SU'  I.M-I'-SIJ'I.H'  >  .<•  '  EU  ■  I » H—  1  ’  — EU  ■:  I .  N  >  > 

1 55  0 

EU  ■  I  •  N- 1  ■  =  SIJ'I.N' 

1 53  0 

IF • SLP.LT . 0. >  STOP  "  NEGATIVE  UNLOADING  MODULUS" 

1570 

!U'I.H-1'=  SLP 

1530 

NM 1  =  NSU'I'-l 

1330 

c^r* 

CONTINUE 

1  *  3  0 

D  TM I N=  0. 3*D  TM I H 

1 3  35 

CP  I NI  T=  SORT'  S  L  <  HI. .  1  >  PO  • 

134  or 

1*5  ML 

INCREMENT  TIME 

1  **  ML 

1370 

PRINT. 

1  *o  M 

PRINT."  ♦♦♦♦♦♦♦♦♦♦" 

1330 

PRINT."  MAXIMUM  TIME  INCREMENT  SHOULD  BE" • DTMIN 

1  7  0  0 

PRINT."  ♦♦♦♦♦♦♦♦♦♦" 

1710 

PRINT, 

17  3  OC 

READ  MAXIMUM  TIME  AND  TIME  INCREMENT 

1  74  0 

READ • 33 , 1 0£ O’  TMAX , DT 

1 735 

DO  £000  1=1. MEL 

173*  £000 

DTD • I =  DT 

0NED3F  Com  I  OS  5*31  uSit'SI  FILE  PftGE  ND.  4 


i  r  r  ii  ntime=i 

1rx1.11;  F’Ehii  feint  find  plot  inteevhl :  hud  no.  of  pt:.  to  be  saved 

17+n  EEhD'  ?  =  .  1  fi£i.f  NPENT.HPLTI.FEEP 

l  T'V  IF  •!  EEP.EC'.  O'  '30  TO  4000 

17*>  PEftt"  lMc-O.  ■NFEEP' 1 1  =  1. KEEP  > 

t •Son  npt;  =  TMH  •  DT+FLOFiT -MELTI  *  i+1 

1  ;3  i  13  n;  ICE=  met  .  ♦e.+F  EEF'+V" 

ism*  if -m:  ice.gt.3Voi*.  :top  "  plot  file  too  big” 

lxl"  4  mm  NFL  T  =  1 

lScO  PPINT.  ••  NO.  OF  MFC  3  P0IMT5  =  -.NMP:  • 

1X3"  FEINT.  NO.  OF  ELEMENT!  =".NEL 

1 8  35  NM I  D=  NEL  3+  ! 

1  :-:4  H  30  CONTINUE 
lSie.ii  DO  34  1  =  1.  NEL 

txro  34  l:  ign-  i  •  =  " 

IX.'S'ii  IF • NPL T. GT . 700 •  '30  TO  135 

1  X'm.n; 

IXOOC  I.OMPNTE  NEI.I  hCCELEEhTION 

1 9  1  oc 

1 3X0  DO  e"  1=1 .NMh: ; 

13*  3"  LN=NLH  .  ■  I  • 

1  '34  0  IF  •  I  .ME.  NMh:  !  •  '30  TO  46 

1950  130  TO  .4*.  40.4V  •  MBTYPE 

1  96.  OC 

i ?7 oc  filed  bottom  eouhdmey 

1  XX  OC 

1  'i'in  4  0  Fd!  c  ■  l .  I  .  =  n . 

.-Om'i  '30  TO  SO 

3  0 1  OC 

3030C  TEHN'M]  TTIN'3  BOTTOM  BDUNDFlPY 

3  o  3  in'. 

30  ■■'*<  45  E  1  =  EF' !  •  c'.  MEL  • 

3033  IF  .HE:  'El  '  .'3T.TDEP!  >  GO  TO  46 

c'O  33  E1=E3 

3  043  FHCT=  1.0 

3  04  3  C’O  TO  4  7 

3044  46  IF-lFLGB.EO.il  GO  TO  47 

3045  DhF'3=  '  !  IG'S'.NEL  '♦!  I '3  •  3»  NEL  '  •♦0.5+iEPi  •£• NEL>  -EPS  <3»NEL>> 

3 046  IF • EP ! • £• NEL • . L  T . EPS < 3«NEL> '  I FL6B= 1 

3  047  FlE3=  hE3*  DhF'3 

3  04:3  F  HC  T  =  FlE3  HF'  1 

3050  47  IiV=  !  NET  >  XM1 ♦FAC  T  EO'  +  'El-Ec'i 

3  06 0  E3=  El 

3150  GO  TO  60 

3 1 60C 

317 OC  FEEE  BOTTOM  BOUNDFlEV 

3 1 S  OC 

3 1  X"  46  FtC C  '  1  •  I  '  =  F  •  1 1  FlMFi 3  3  '  I > 

3300  60  CONTINUE 

33 1 0  I F ■ MOD • NT  I ME . NPENT i . ME . 0  >  GO  TO  101 

3 3 c'O  FEINT.  '1  TIME® '.THYME 

33  30  PEINT.  1 16-  1*  l.J. EPS-  1  •  1  -  *  SIG  <  1 «  NM1I"  »  EPS  •  1  >  NMID'  .  SIC  ■  1 .  NEL)  « 

3331  ■: EPS • 1 .MEL' 

3340  101  I F ■ MOD ' NT  I ME . HPLT I * . HE . 0>  GO  TO  103 
3345  IF -FEEP.EO. 0-  GO  TO  103 


une 


£  £  r  o 


c,c,74 
££*11 
£  .*• » > »’» 

£  j  l  u 

£  3  0 
£j4».'C 
£  5**  ML 
£  -:»r.ML 

£  5  r  >!• 

£  *:y  0 
£  5'40 
>4Sn 
£4*0 

£4  roc 

£4&  ml 
£4'r*  Ml* 
£*.00 
£5  y£ 

-,ei  i'iS 
£*•  1  0 
£M5 
£*5 1 7 


£*,*=»() 
£*00 
£*  1  oc 

£*£0C 


7 


,  []ht  1 1 1 :  S:  VI  OS  £U  SI  FILE  PHUfc  MG. 


TIME 'HFLT • =  THYME 
DO  1  "S  1  =  1.1  EEP 
mi=  'i-i  '♦hft:+  mf-lt 

IIS=  •  I-l+> EEF'*HFT:+  MPLT 

n:--=  ■  i  - 1 +c*f  eef- •  ♦mpt  :  +  mf-lt 
H4  =  >  I  -1  EEF' '  ♦HPT  :  +  mf-lt 
Mc.=  ■  I -1  +  4»1  EEF  •  ♦MF'T  :  +  MF-LT 
I  =  Mr  EEF  ■  I  ■ 

FF-LOT  1  H  l  ■  =  1 1  i?  '  1 « 1  ■ 

FPLOT'MS'=EF-'  '1.1  ' 

ff-ldt  '  m;--  ■  =  rt. c  'l.r  ' 

FF'LOT  •  H4  '  =VEL  '1.1  ' 

1  0£  F F'L Cl T  ■  r,f.  ■  =DI  : F'  •  1 . 1  ' 

MF'LT=  mf-lt  +  i 
in}:  M  T  I  ME  =  MTIME  +  1 

do  4?  i=i.mmr:: 

4  -i  F  '  I  '  =  0. 

THYME =  THYME+I'T 

1HTEUPHTE  hCCELEPhTIGM 

DD  5F.  1  =  1,  NMFt ; ; 

VEL  '  1 .  I  '  =  VEL  '  c .  I  ■  ♦  HC  C  '  1  •  I  •  ♦I'T 

n  i :  f  ■  l .  i  ■  =  D  i :  f-  •  s .  i  ■  +  vel  •  1  •  I  >  ♦M 

SF.  CQMTIMIjE 

I F  '  MET-,  PE .  EC' .  3  ■  VEL  •  1  « NMfi 5 3  '  =  VEL  .  £ .  MMH3  3  •  +I'V 

COMPUTE  IT PHI  M3  HMD  3TFE33E3  > 

DO  1  Of.  1  =  1.  MEL 

ru:nT=ij 

ii  =  i 

LM=  MLh','  >  I  • 

HYTHH*  ETH ■ LM • 

DEP }  =  D T ♦  •  VEL  '  1 .  I  '  -VEL  •  1  *  I  + 1  >  >  '  2  • I  '  -Z  •  I  + 1  >  • 

I F  '  HE }  ■  DEP  }  '  .  >31  .  TDEP  .  '  '30  TO  58 
I F  •  HE  :  ■  EP I  ■  £ «  I  •  •  .  '3T .  T DEP '3  ■  '30  TO  58 
VEL' 1" I '=  VEL • £• I ' 

IF  •  I .  EC'.MEL  '  VEL  '1.1  +  1'=  VEL'£.I  +  1' 

DEF  :  =  '.'. 

5S  EP;  '  1  •  I  •  =  EPI  •  £*  I  '  -DEF'i 

I F  •  EP  :  ■  1  .  I  '  .  LT  .  EMH: : '  I  ■  '  L  3  I '3M  '  I  '  =  1 
EF=  EPI  .1.1' 
e:=  EP;  ■£.  I  ' 

IF • EP;  ■ 1 .  I ' . 6T.EMHS • I > . HMD. E3 . LT. EMHX ' I >  >  L 3 I6N < :  I >  =  I 

E . M 1 =  ep: 1 ?« I '  1 

::=  :ii3'£.i' 

: :  m  l  =  :  1 13  •  3  •  i  •  | 

dth=  o .  ' 

DTM=  I'T  | 

E  I  =  EF 

I F  '  L  3 1  '3M  '  I  '  .  EQ .  0 >  '30  TO  6d 
IF  •  EP  3  •:  1 .  I  '  .  '3T .  E  3  •  '30  TO  4  00 

UHL OH D I MU 


■  -  .'••■a*.  >- 


LINED. :P 


rOMT 


PILE  PAGE  HO 


1  0: 


?1 


Oh  ct-  31 


aa mo 


£t~.  3 5 

mmh::=  roj-LM'-i 

C*t*  38 

IF  ■ES.LT.EINl  ■  GO  TO  300 

£6  4  0 

I F  •  E  : .  L T .  EHAX  ■  I  •  ■  GO  TO  £35 

ihSii 

ENE  •  I  •  i  •  —  emax  •  I  ■ 

£if.*0 

no  £35  j=i.nmax 

£6.  r  o 

;.OC 

IF  •  IMAX'  I  ■  .GT.EU'LH.  .1'  ■  GO  TO  £37 

at*  8  0 

£87 

dummy®  smax  •  I • 

£*85 

I  FLAG ' I '=J 

£680 

jma:  :=  nma: j+£ 

£700 

DO  £30  HHE=  c'.JMAX 

^705 

div=  :u 1 LN«  J • 

£  r  06 

I F  •  AM£  •  LH  ■  .  L  T .  1 .  EG  ■  D I  V=  D I V+  AM£  •  LN  > 

£710 

ENE  '  I  *  HUE  1  =  EHE  •  I  .  HUE  - 1  —  '  DUMMY— EU  '  LH  ■  J  ■  •  HIV 
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: :  m  i  = : ' 

£  78  0 

3  0  O 

: :m i  =  Hfti  •  ii.es. ef.ei.nmax,  i.» 
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F'EFlD  1  3  'r* .  X  0 i1  -  ii.  i'ii.  MnL E  *  _ FX.  IF  Y * ; CL .  YL 

d:  ;=  :f:  : 

D't=  :f',' 

Ih=  '  MbLE  +  l  •  ?rj+l 

THETR=  Ni5LE*3.  141^6?.  18m. 

C  =  CO:  I THETFc 
:=; IM'THETR. 

Chll  plot  ‘  :»;o.  y'm*  - • 

IF  i  I R .  ME . £ . HMD.  I Fl.  ME .  4 '  MO  TO  E 


b 1  0 

I DUM=  1  LX 

6  c.  0 

il::=ilv 

6  3 1.» 

ILY'=IDUM 

rr-4  0 

::dum=  d:: 

if.  5  l.l 

D:  :=  I'Y 

6k-  0 

Di'=  :DUM 

6  7  0 

::dlim=  :-:l 

hXl) 

XL  =  i'L 

K.’-*  0 

YL=XDOM 

7  0  i’i 

E  MC'/.*  Ml  •  I  LX 

♦  1MC:  .  >  IH 

710 

Ml  ('=  Ml  •  I  Li' 

♦  I  ML  ‘V  '  I H 

7£  n 

:v::=  d:>i:v: 

■  I H - ♦XL 

7  ?  0 

rvv*  DY*i:y 

•  I H  ' ♦YL 

74  0 

D:  :=  d:>iix 

IH' 

750 

DY=  PYMIY 

IH. 

7  if.  n 

x.  :  l ♦re: : 

IH  ■ 

770 

;*  (  i  -  vl*h;y 

IH  ■ 

7  7  c 

:  C=  XT  ■  in •  ♦: 

L 

QHk  D  CP 


CDNT 


1  i.i : 


VI 


FILE  F'HHE  HU. 


x  ii  n 
1  I' 

x  £  ij 
1 1 


*4  M 

_  c  .  . 


1  Ii 


X  ,  II 


.  T=  i  T  •  IH  ■  ♦YL 

chll  : ,MFOL • : z 
r  mll  m:  : i  :  i  3  •  z: 
chll  h:  i :  l  ?•  o. 

ChLL  FLOT  ■  ii.  •  ii. 
DO  Hi  1=1 «M 

n=  .  .  i .  *c  :f:  > 

v=  ■  l  :f:> 

C HLL  FLQT.H.V.£ 
com  I  HUE 
CHLL  PLOT 
PETUPH 


_  ,  .  0 .  .TITL.hT.  IH' 


.  0.  «  LHfL  •  il;  : 
~YY  *  LHfL  *  I L i 


•  no 

.  ml 


.40' 

.  os  • :  l  < 

.  ii x  •  YL  ■ 


mop  ■ 

MIip  ' 


il:: 

1  Li‘ 


1  .  •  1  ■ 
1  .  .  1  ■ 


:fy' 


I 


li.  -  VO,  -  i 


:X>;  i.i 

m-xiiJ 
I.I  0  I.I  1 
01  01 
01*  01 
II  £  I  f) 

OS  ll'j 
l.lix  Ii') 


£00  FOFMHT  ■ 

EMD 

s  i J ;  E  S  .  i?  T  L  I  T 
:  EXECUTE 
:  l  i m i  t  : :  :-.o. : 

:  F  FtTFL :  iV .  U .  L  •  P  0 1 D  JOC  FXDIiHT H 

:  THF-ET:  0?. X£I . ♦CPL 

:  IMCOI'E:  I  I  MF 
:  EMD  .'Of 


'01 
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'  K  1  . 

ci  •  >  *,  .. 


L L 


0NED3PD 


4.  0NED3PD  is  the  driver  program  mentioned  in  Part  V,  which  can 
be  used  in  a  tr ial-and-error  process  to  establish  values  of  the  viscous 
parameter  n  and  the  spring  .  This  code  is  presently  set  up  to  run 
on  the  WES  DPS/1  time-sharing  computer  system  and  makes  use  of  a  data 
file  which  is  very  similar  to  that  required  for  0NED3P;  however,  because 
the  time-sharing  system  is  interactive,  many  of  the  input  quantities 
required  for  running  0NED3PD  are  input  from  a  teletype  keyboard  by  the 
user . 

5.  Referring  to  the  data  file  shown  on  Figure  A1  for  0NED3P 
batch  runs,  all  that  is  required  for  a  0NED3PD  data  file  is  the  problem 
title  from  Section  1,  the  static  stress-strain  curve  (M^  function)  from 
Section  2,  and  the  forcing  function  from  Section  3.  All  other  input 
information  is  asked  for  by  0NED3PD  at  the  time  of  execution  through  the 
process  of  interactive  questions  and  answers.  A  typical  0NED3PD  data 
file  follows: 


"  r'HJ  IM. 


0  0 

7  t.-  . 

.  c 

: 

c  :■ 

.Hi 

< 

Cj  ; 

.  1 1 

o.  . 

1  0  5 

J 

■■  LOhHi 

FUHC 

non 

MF  h-m : EC 

cl 

0. 

0. 

cl  0  5 

0 .041 

0.474 

0 . 

0&t' 

0 . 7  54 

0.1  os  / 

1 . 

T  *5 

1  50 

C. 

07  4 

0.  l='i 

2.  177 

0. 

1  74 

4 .  04  5 

0.  177  1 

1  5. 

'■* 

7  1  r* 

17. 

5  M 1 

0.  74  1 

1’?.  54}: 

0. 

7v. }: 

£0.641 

0.75*:. 

7  7 . 

5  ^7 

5  5 . 

5r>  7 

0.  2-  0 

cr 

0. 

557 

0.  574 

47. 

:  r. 

1.1 . 

•:  *7 

4  ? . 

1  :-4 

u .  4  l  f* 

40. 75’? 

o. 

44  1 

i  r .  ■?  i 

0.  45  5 

751 

4  -  »:■ 

7  j:. 

1 

0 . 5  04 

70.  754 

0. 

5  5 1.» 

14.  ?$? 

0.  557 

Id. 

54 

i  j , 

574 

f }  r.  j: 

c  . 

757. 

h  .5-7 

174 

0. 

417 

c  — 

0  .*=41  J 

}  TITLE 
STATIC 

CURVE  STRESS-STRAIN 


LOADING  FUNCTION 


6.  A  listing  of  0NED3PD  follows. 
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■  fit  ■ 


•i  ['  >  ! 


1  ii :  It.:  t *3 


04--20  ! 


FILE  PARE  l!C. 


1 


■■  1  1  i  It  1-44  I  FLOT  -  •  F 

'T  'I" 

ir,':-  ClMFriSF’D —  HP  I\-ER  .-nF  GNEI'VP 

'  4  mi 

it.  niiE-I.lHFti:  I  OUR.  EEHRVIOF'  QP  Fi  POINT  OF  VI .L GELA.TIC  MhT  4 

THE  '"  Di  v  T I  Tl.'T  I vs  IEHhVICP  MRV  BE  PEFFE  SENTED  IV  R  MECPPr-I 
I'. MODEL  LONS  1  STINE-  DF  h  PRFhLLEL  SPFIM3  HUE  MAI  "‘ELL  ELEMENT 


•  ii.  .*’•*:!*  UNITS  RYE  ANY  set  DF  COM!  I  I  TENT  UNITS;  PDF'  INI  TRUCE, 

ill  i.ti.tr  MEHFiF'H 3 ORE 3  •  MS  EC»  METEF’3,  Kt» 

1 1 1  l  fnl 

iilc'i  CHRP  AC TEF ♦ 1  AN  1  *  T  I  TLE^F-F-,  F I  LE*S4.  F'T  I  TL*40 

ii  i  *iii  C QMMOM  1L  *50  •  *  I- U  >' 50>  »  EL  •  50  >  •  ELI  <50*  *  ENE  <i  0>  •  IFLA<5*N."L* NSU 

ft  i  I  fi  0 1  MEN  Fieri  TIME  <  5  0  0  •  .  I.IJM  1-50  0  :•  ,  IUJME  <5  0  0  >  .  EPS  •:  3  >  *  S  I '?  •  3  > 

M  J  7*.  Ill  MENTION  FFLDT  •  1  0  Fi  1 1 0  > 

l;  !':'*  !  l.HLL  F  PRr'AM  ■  1  «  SH  > 

i'  ■  * *’•  i_  fill  f ,\df t  <6  7 ,  i  *  i  »  o * 

Cv  Mm" 

OS'  1  fi  PRINT. 

ii-'' S 0  5  PRINT,"  INPUT  FILE  NAME" 

lie?  ft  PERU.  FILE 

o.34i.i  EMC  DIE  '  FILE*  1000* 

nS50  CALL  ATTACH  '  3 :j*  F  IL.E>  3*  «  » .> 

Opt.  0 
03. 7  0 

Mi 3  0  READ  •  33.1  OS 0  >  p T I  TL 

03  3  T  PRINT. 

i.i  .-  I'm  F*  F‘  I N  T . 

Mill'.  PRINT. PTITL 

0 ' 3  •  PRINT, 

-1  DO  10  1  =  1*3 

0:70  SIi?'I,=0. 

o ; 3  t  eps • i *  =  o. 

''.'"r '.  10  CCNTINUE 

os sma:  =o. 
i'.c  i  o  sm=  o. 

•i/.:  o  emr;'=o. 


31=  0. 

!  «S4  »'i 

3S=  0. 

ii5e*0C 

.17*7’  0 

PRINT.  ••  ETA.f 

S" 

! 1  .7*  •  i 

PEAR  ETA.  RMS 

i'  f.  r- 

F'E  AO '  1  Oc'l." 

NSL.N'U 

fi.r  r*  i.i 

PERI' '  3  3*  1  03' 0’ 

■Su  •"  I  •  .  EL  fi  >  ,  1  =  1  .NSL.i 

J  ■  » 

F'ERTi  *  33,  1  03  O  ' 

• 3U • I • *EU • I> • 1  =  1 • N : U  > 

1 

I'UM  1  •  1  •  =  ii  . 

1  : 

RIMS  •  1  -  =  0. 

1  ;4.» 

NC'LT=3 

1  **; 

i  •>'  iv: 

IN1I.IT  SUPFACE  :T='E33  time  histgpv 

1  .*•  1  * 

1  »' 

THYME.  =  ft. 0 

1 

:  1 1?  ■  i  •  =-  sstps 

'  I  H  i' ME .  1  * 

H- -’V 

14  *•.; 

ESTRtL 1 : H  R 

TA5LE  OF  MODULI  FDP  THE 

FUN'S  TIC- 


<J)  u  I 


•  r  OIPO 


.[IM 


1  0 :  It.:  1 1 


04  £lj-?M 


FILE  Pm-.E  NC 


I  -*  7  •  *C 
1  4  ;•  •  - 
1  4 «■' 
14.M 
14 

I  4  '** ' ' 
i  S M 
1  r-»4  !  I 

! *■  **•  i  t 
i  c.*  l'i 

i  *5  r  0 
r;  >  »:■ 

172  o 
174" 
j  7*? 
i  7  *  0 

\  7*0 

1770 
1  7  :-M 
1 0 
i M  l : 

1  r.  1.1  c: 

J  V'->* 

2  ?•'»»> 
1*4  0 
J  V7ri 
1  ?  0 

1  S OC 


:vc'4 


c’cV^O 
I*  2: 0  0 

*  •  2:  r» 

*  T'4  OC 
74  7  OC 
74  ;.  0»2 

£’  4  *  or 
Jf.  07 
j*?'  o* 

-  •.  i  =. 

CY  T.C 

£'.4p 


j>,  or. 


on  £4  ii=£.n:l 

'  l  p-  •  •  M  1  -  '  L  M-  1  ■  •  -  EL  *  N ‘ -EL  *  N— 1  -  ' 

IF-  -LF.lt.  0.  ■  ITOP  "  NEGATIVE  LDrtMrc;  MOOULU:  " 

£4  .  L  '  N-  1  •  =  :  lF 
n;u  =  mil  -i 
OG  EG  M=c  •  rl :  u 

.LF'=  ■ ilj'li-l --:u-h-  •  -EU‘N-1  --EIJ-N- > 

f'.i-M-l  •  =  '•  i1 '  f • ' 

IF  .  CLP.LT.  0.  -  *  TQF  "  NEGATIVE  Ut'n_C*PING  MOBULU • " 
EG  CLI-N-I  •=  -LF' 

'i  ei  i  =  h;l<  -i 

PRINT."  Mm  :  I  MUM  TIME" 
remd*  tmm: 

F  ~  I  NT  . "  TIME  IMFFEMFMT" 

FERC.PT 
OTG  =  Of 
NT  I  ME  =  1 

FRINT."  FEINT  INTEPVhL  F-NO  PLOT  II-TERVE-L" 

PEFtB«  NCPNT ,  MELT  I 

NPTS=  THmY  ■  01*flChT  .NRLTI  •  -+1 

n:ize=  hrt;*£+^oo 

IF-N:  IZE.l’T.  lOOf-T  •  IT  DP  "  PLOT  FILE  TCO  BIG" 

NPL  T  =  1 
■30  CONTINUE 
L:IGM=  0 

I F  • NPL  T .  GT .  1  Ti  0  0  .  -.'0  T  C  1  Elf. 

IF  •  MOP  -  NT  IHE « Nr’PNT  •  .HE.  O'  GO  TD  101 
PF’  I  NT  .  THYME .  I  1 6  -  £  -  >  EP  i  ‘  £  • 

101  I  r  -  MOP  NT  I  ME » NPL T I  >  .  NE .  0  -  GO  TO  103 
T  I  ME  •  NF'L  T>  =  TH  YME 
Nl=  £♦  NiF'LT-  1  -  +1 
N£=  1 1 1  *  1 

FPLOT •  N1 - =  IIGy£> 

FPLOT  ■:((£.:-=  EFM  <£> 

NC'L  T  =  NPL  T  +  1 
103  NT IME=  NTIME+1 
THYME =  THYME +0T 


COMPUTE  NEW  STPMlN 


;<0  f>  T  =  0 

:-■=  :  i g - 1  • 
h  YTF-N=  ETF- 

IF-EPI -p -.LT.EMPY-  LIIGN=  1 

e:=  ep: <£> 
e  s m i *  ep : -  3 • 

i  f  •  £  i .  G  r .  emh:  .  hno  .e3mi.lt.  emh:  '-  l  :  i  gn=  t 

::=  :iG-£- 
: :  n  l  =  :  i  g  •  3  - 

i-i  n=  r-T 

IF  -  L  .  I  GN  .  E1.'.  1  *  MO  1  0  c-ir. 

IF  - E ’ . GT . £ i M 1  •  GO  TO  400 
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,  f  -*  A* 


- ? "rnmtmmmmmmmmmrn 

(.  »i*S! *  atrr*. 


1  (i :  ! F. :  1*  II 4  £ it  «c i 


FILE  PhfE  HD. 


'i.iLD-P  Il>6 


:iMh:.=  fi'u  -i 

i f  •  e : .  l  i .  E : m .  hi<[i .  t -  fn . L t . emf*:  •  go  tq  ?.*>■ 

E>-F  .  i  .  =  EMh' 

LiQ  iih'f  J=1.('wh; 

IE*  .  ,-,I H ,  .  I.. T  .  £1.1 1  l  ■  •  lif'  TO  Go 7 

iL-rif",  =  :mfx 
1  FLhG  =j 

.*  H.TGX-J+G 
PD  r.' 1  r* 1 1  >v‘*c.=  c  *  -_*MF>N 
P  I  V=  .  I.'  *  J  * 

IF  •h;',c'  . LT.l.E*'  IHV=  DIV  +  Fii-1£ 

EHE  •  Mi*c  '•  =  EMF.  aiHE-l  •-•I'l.M'-iV-Fl.i.  J  • 

PUMMY=  £IJ  •  ■ 

.>=J+1 

"M  l  =  (ir- !  ■  £  I  M  1  ,  £  S .  NMP  V  %  1  > 

EF  =  El  PI  *XM1  .  pME.HYTF*h.:3SNl  •  .  I  ■  '  r  .  E  EM  1  •  E  S  •  pTH  •  IiTD •  P IM1 > 
:  m  i  =  :  [•’  i  +  n :  m  i 
GO  TCI  70 

PELOFtp  IMG 


■  i  »:• 

4  0  0 

CDnT  1  HUE 

-  •  •-  o 

HMh'<=  MIJ-l 

-.0 

4  1  0 

>:*.(  1  =  PM !  •  E  :  M 1 « E  S  J  ftMH  •  E> 

l'|  7*  f  j 

4  If- 

£F=  E  IP  1  <XM1 «  PMc- •  ftVTT*M«  J ■ 

:mi=  emu-  p:m 

I.-4  !  f 

4  c.  0 

IF.iEF.LT.EMP,:>  60  TO  7  0 

}:41o 

EMP  =  EF 

*:41  7 

INh  =  ?F 

34  £0 

uC  Tf)  70 

.:44  •>.: 

VIRGIN  LDhPIMG 

54*  »*.£ 

3:4<-  0 

C'C 

CQNTIMHF 

HMFi.<  =  H.L-1 

>?  i  0 

o  j. 

XM 1  =  HM 1  .;E  :  M 1 .  EE  >  MMHX . 

f*  >'• 

64 

EF=  E I F' 1  *  XM1  *  RMS'*  MVTftHi  3  E 

:•>.  • . 

:  !•>  i  =  :  hi  +  tsMi 

.V  l  o 

IF  >:£F . LT. EMPXi  GO  TO  70 

j 

of 

:>i  hx  =  -t 

:-  .r  "4 

EMP/  =  EF 

y  .  ^  ;  r 

70 

fOHTIMl.it 

4  1  > 

'  I G  ’  3  •  =  ” 

4  )  K 

.  I G 1  £  1  =  .SF 

4 1 1 

EPS • 3>  =  ti 

4J 

£F- ?*.£•>=  EF 

4  1  .  ‘ ■ 

PTG  =  DTm 

M  4<  ’ 

4  )  r  ?v 

j  _  , 

h>'P  .  U.'-FPi'E  ;  TP£  ;  ; 

' 

: i g • i * =  : : tp ; . th .me  >  nr » ME 

i  1-  • 

IF  .HPIT.LT.  lOOv.PI-p.Th.  ME 

4,. 


*  *•  -*5  A  Mr  jq  £  ]f  * 


4 1. 

niiiu  nutft'i 

■>- 

PLOTTER  GUTP"T 

.i,-  " 1 

i  .  *< 

ice. 

Ifi  M  .  ‘  !<:;  i'Qi.1  WANT  F  L  D T  T F  P  OUTPUT t " 

4,.  '•< 

f>  a;  . h.- . 

\  -  1. i 

I  r  '  pj  -  .  Mr  .  "  V  "  •  MG  TG  It  11 

4  1  »* 

i  *  ;*> 

r«4'LT=  HPlT-J 

f  p i mt .  ’  next  plot - 'vpe  cope" 

•♦*4 

reap.  n  i 

;  -*  .* 

I F  •  1*1  .  cO.  o  •  '?□  on  160 

4*1  -*■. 

60  TO  •  1  3E*  1  3S«  1  33  •  .Ml 

4  ■  *  < > 

I  ;*J 

r-,c  i  r .  '  •;  TRF  3 3  -T I  ME  " 

4-  1  0 

do  t :•  i  =  i<hplt 

-  ;  | 

ii=  l*  •  1-1  •  ♦! 

4  '■  •••■ 

1  *•'  -■ 

Kii'M-l  '=  Fr  LOT  <14  • 

CALL.  FLDTT  ■  TIME  •  1<  *  I'UMl  •  1  •  .MELT.  1  »E«PT1TL 

4  ft 

GG  TCI  1  '.0 

1-  T  •» 

1  3S 

PRINT , "  .TRAIN- TIKE" 

4  ;  "  11 

Tin  13-  1  =  1*  MEL  T 

4  "■  ;r  0 

r<=  £♦ > i - l • +E 

4f.  ^  ft 

1 

rn.iMl  .  I  •  =  FPLGT  <M  •♦100. 

44.  i  i*i 

CALL  PlCtT  •  TIME.  •  1  >  .  PUMl  •  1>  •  MPLT.  1 . 3  <  P  T  1 T  L 

4-  ft  * ' 

AC  to  130 

4r.  *•  ,  : 

1  2:  5 

PP1MT. "  3  TRESS-STRAIN" 

4  T'  l  1  1 

I'D  13?  I  =  1  <  MPL T 

4rc.fi 

H=  .?'♦•  1-1  >  +  t 

4r-" 

hn=  n<T 

4-  ro 

DOME:  •  1 1  =  FPLCrOM 

4r  ft 

i  :-■? 

PUM 1  • :  P  =  F PLOT  'TIM  ■  ♦  1  0  0 . 

4  ft  i> 

call  F'LQtt  <pumi  < i  >  *  dune •  1  >  *  mplt •  ?.  e«  ft  i  tl 

4-1  1  0 

GO  TO  130 

4  4  “•  r. 

160 

PEW  I  HP  3'? 

4>L*  4  ft 

CALL  PET  AC  HO'?.  .  > 

T.t'l  4  ft 

PP1MT . "  DONS-DONE-PONE" 

ft  ft  ft 

PRINT."  pa  YOU  WANT  AMO t HER  RUN7" 

r.  0  4  ft 

REAP.  AM 3 

Aft?,  ft 

IF  CAM-:  .EO.  ’  Y“>  GO  TO  S 

A  ft  '  I 

STOP 

r  ft  7  ft 

1  U'iO 

FORMAT CTE4. " • "> 

A  ft  ?'■  1 1 

i  o 

FORMAT • V> 

ft  ft  '*4  ''I 

1  0  ;•  >3 

FORMAT  1  I  1  OOF  1 E.  S'. 

ft  1  ft  'I 

!  04  0 

FORMAT  >  ••  LAVER".  J5>  "  V  ICC  07  ITY="  .FIS.  3.  " 

-•1  1  ft 

1  O'-'!' 

tf OP r'E’T  1  7X.  "MOPE  "  •  6X»  I G"  • ’«  "EPS  "  »  9\t  "AC 

ft- 1  ft  •:> 

t  l'-  l 

FORMAT  •  1  1  (i *  4F 1 E .  S •  F 1 E .  7  ‘ 

-  i  2- 

EMP 

-  1  4  *V 
*1?»* 

FUNCTION  3iTRTOT.rO 

ft  *  r_>  " 

CHARACTER  T  1  TLE*r'r 

ft  i  r  ft 

PI  MEN:  I  an  >.‘600.  .  Y'.E.OO' 

*  i:f  •• 

IF'M.GT.  1  >  GO  IfJ  10 

••  •  J  ’4ft 

REAP • 3?« EOO .  TITLE 

r  ft  ft-  ft 

RE  AO  i  3  ?  >  Eu  <•  *  MPTC 

ft  ft:  i 

Pc  AP  •  3  ?  •  E  0  0  >  '  t'  ■  I  •*>.•!>»  1  =  1*1  IF’?  7  > 
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COM 


1  (isle-:  It  04  0-81 


FILE  Fm.SE 


1  4  i;*r 

1  4  ?.  lj 

I'D  34  i)  =  £.N;L 

1  4  f ! 

: lF'”  ■  ;l  •!" -  l >n-i  ■  ■ 

1  EL  1  N  '  -EL  ‘  N- 1  >  :• 

1  4  ’  II 

if . :lp. lt . o. ■  : top  " 

NErhT I VE 

LORI. IMG  MOTH.ILU I  " 

14.ii 

34  . L  <N— 1 • s  ILF 

1  4  ‘4 '» 

nmi=  n;l  -i 

1  *:•  »'i 

tin  3'F-  N=3  >  N  .  U 

IMu 

ILF-  1  1 1.1 '  It  - 1  '  -  ill  *'N  ■  > 

'iEU'N-1 ■  ■ 

-Efi  ■  N  '  > 

!  *■ i  t 

fn  • H- 1 •  =  : U ■ N ■ 

lFiilP.LT.Ci.i  'TOE'  " 

NEGRT IVE 

UMLCRPING  MO PUL U 

\  sro 

.-if.  :u'N-i'=  :tp 

1  0 

HM1=  N : U  -1 

i  r :  »> 

PRINT. ”  MR  .  I  MUM  TIME' 

1  .“4  0 

E'Ehti.  TMr*V 

17? '1 
'  7f.  0 

1  7f.f. 

1  770 
1  7  3  0 
179  0 
1  9  0  0 
1  -0  Me: 
ICO* 

1  3  0  0 

i  7:4  M 
IV  70 
1  :  SO 
!  ?  0  Of 
33  1  0 
33'3'  0 
3  3  40 
33'?0 
:v‘?4 


3 j  0 
'  4  Of 
”  4  7  Of 
14:-.  Of 
14  9  Of 

■Si:.-; 

■SI'S 

=.l? 


3  0 


FEINT. "  TIME  INCEEMFNT" 

EERP*  I'T 
DTD=  HI 
N  r  I  ME  =  1 

FEINT. "  FEINT  IHTEPVRL  RHP  PLGT  INTEPVRL" 

EERP«  NCPNT .  MELT  I 

NP’T'I  =  THRX  ft!  ♦FLGRT  (HPLTI)  >  +  l 

N‘.IGE=  NET  :  ♦E'  +  SOP 

IF'M:  I2E.RT  .  1  000.-  •  .STOP  ••  PLOT  FILE  TOD  PIG" 

HPLT-1 

CONTINUE 

LI I UN=  0 

IF  .'MPLT  .  GT .  1  000'  CD  TC  IE'? 


j  :•  l.< 

'r-  0 1 . 


1 F • MOD (NT  I ME  .  MEPNT > . NE .  O'  GO  TO  101 
PEI  NT  .  THYME  .  -  I G  ’  c.’  .  EF‘  S  <c’> 

101  I F 1  MOP ■  N T I ME . NPLT I > . HE . 0>  GO  TO  103 
TIME  'NF'LT>  =  THYME 
Nl=  3*  UlELT- 1 1 + 1 
H3'=  1 1 1  + 1 

pF'LOT  '.Ml  '  =  S I G SET 
FPLOT  <N-2>=  EP-?  <2> 

HDLT=  NELT+1 
103  NT IME=  NTIME+1 
THYME  =  THVME+PT 


COMPUTE  MEM  -ITPRIM 

NC  M  T  =  0 
■e=  :  JR.. i  ■ 

R.TRH-  ETR 

if  .EE : ■ c ■ .lt  .  emp:  '  ■  l:ign=  i 
e  :  =  ee  ;  •  iv 
e:mi=  ee:-?.. 

1 F  <£f . G I . EMh:  . RNP. EIMt . LT. EMR>"'  LI IGH=1 
:iG'E. 

::mi=  :ig<3. 

[iT  M=  T'T 

ip'l:ign  .eu.c.  so  to  re 

IF  'E ; . GT. E I M 1  .  CO  TO  400 
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OCNT 


1  0:  18:  It- 


FILE  R"H*5E  NO. 


04 -  £ 0  SI 


' ;  o  D£®  •  _  IP1-;  I  ■  *fi 

.  t  •  *>  d:<=  <.:i-simi>*b 

■  V  MI.  D4  =  •  £  i  — e  i  m  i  • 

’.-•10  EIP1  =  EI+D£-D?+D4 

; Vc.-  0  p  S  K !  =  C  1  ♦  •  E I F 1  -E I  ' 

rET'.O  RETURN 

EE 40  END 

7E  '.or 


7  c' 8  00. 

"7£ E 0  Sl'PPGUTINE  PLOTT  *.X»  V»H»  ILXT.  ILYT»TIU.> 

"738.0  CHfiRfiCTER"  *1  HN 3 , Lfif L* 1 5 

7£?0  DIMENSION  X<  1>  ,YU)  *LAtiL  <7>  .NO  •..?>  ,HI*P  <?">  »N2X  <'4>  «NZY<4) 

7  7  DO*  DIMENSION  INCX<4*  .  INCY<4» ,  ISV>.  <4»  •  ISVY<4>  »  I IX  <4  >  *  I  I  V  <4 

7  3**5  DIMENSION  fiT  •  4  >  ,  XT  •  4  .  ,  VT  <4> 

' 1  *.*  DhT  H  !  N'.  K/-l'  *N  1 »  1  /»  I NCY  1»— I*- 1'lc.I  S  ’’y. y  0  •  1  ,  1  »  (1/  9 

.’?.i  "■  isvYvo»  0- 1  •  ix,  i i:<x  1  •  -1  -  —  1 ,  i/»  iiY/ii  1  , -1  , -\* 

r  3 3  0  I’R  Tfi  NZM  ■  0  ■  ,  r  18Y  ■  0  ,  0*,  —  1  ,  —  1  / 

7  3 DR* Tfi  H  T  'li .  ■  “•  U .  »  li .  .  E  7 1) .  'XT'.  7  ,  —  — 1.3,  .  S  '  ,  YT^  .  "9,  .  7  ,  — . 

' 94 0  DfiTfi  N*.  -  1 1 , 1£,  1  0,9«  1 3  ■  9 «  1  0  ,NDF  '1  *  1  •  1 9  1 ,  1  ,  1  ,  £-‘ 

r.'SO  LHlL<:i'=  "TIME  -  msec" 

738  0  L81BL  <£">  *  "STRESS  -  MP*“ 

737**  LflSLO'*  "STRAIN 

74 Af*  PRINT ,  "  WANT  PLOT?" 

74  v 0  READ, 

7.140  IFtfiNS.NE."Y~>  RETURN 

74:.0  I  LX®  ILXT 

74  R.  0  ILY®  ILYT 

74, -0  PRINT,"  REfiD  XO,  VC*,  THETA,  SFX,SFY,  XAXISL, YAXISL  " 

■  4 0  REfiD,  XO,  YO,  MULE,  "2FX,  SFY ,  XL,  YL 

749  0  DX®  SFX 

7500  DY=  SFY 

7510  Ifi=  <N6LE+1:>x90+1 

75. -0  THETfl=  N*3LE*?.14159£65  180. 

75  30  C=  COS  •  THETfi :• 

7540  3=SIN  (THETfi' 

7550  PRINT,"  DC  YOU  MANY  THE  AXES  DRAWN?" 

7580  REfiD,  fiNS 

7570  CfiLL  PLOTS '""A") 

75  r  0  *  ft  L  L  PLOT  < X 0 «  Y 0 ,  —39 

7590  IF  <ANS  .  HE.  "  Y">  *50  TO  5 

78  00  I F  <  I  ft .  HE .  8 .  AND .  I ft . NE . 4  >  *30  TO  £ 

7610  IDUr1=  I  LX 

78 £  0  I  LX- JL V 

7620  I LY= I DUM 

7t  4  0  XDUM=  DX 

78  5C  DX*  DY 

78ir.li  Dr"®  XDUM 

787  0  xdum=  XL 

7830  XL  =  YL 

78  90  VL=:  "DI.IM 

7700  £  NCX=  NC  '■  I  LX)  ♦  I NC X  <Ift> 

77  )  *:•  r»CY*  NO  '  IL  Y  >  ♦  I  NT.  Y  <  I  ft  ■ 

77 £0  :v  ;=  D;  ♦  I SVX  '  Ifi  '  •XL 

7*  30  IVY®  DY*  I  5  VY  '  Ifi  •  ♦YL 

774..  DX®  DX*  1  IX  ‘  Ifi  • 
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1  n:  It.:  it. 


04  c’  0  •  o  1 


FILE  PACE  MCI 


l«V=  I  I  <  I A  ' 

zxx=  :  l*mz:-:«  ia- 

ZYY*  YL*MZYvIA' 

::z=  ,\‘T  •  i  A'  ♦XL 

VZ=  Y  T  •  I  h  '  ♦YL 

call  '  iHfol  yz- .  or*  t  i tl •  at  < i fl >  .ar* 

CALL  AY  I  :  13‘ZX.»..  0.  *  LAEL  < IL; . «  .NCX» .  1  •  XL  •  nl'F  ULXJ  *  0« 
i.  ALL  A  II-''  0 .  *  ZYY  *  LA  PL  1  I LY  1  *  HC  Y *  .  1  ■>  Yl  *  MI'P 1  ILY  ‘  *  1  * 
5  DO  10  I  =  1  , N 

H-  I >  *c  '  SFX-  Y  •:  I  >  *S'  ?PY> 

V=  <>:  >;I  ■  ♦:  SF»Y  •:  1  1  +C'  iFY> 

CALL  PLOT CH, V.£> 

10  CONTINUE 

CALL  PLOT  <-X0j -YO. 

C  Ai.L  PLOT  C  0 .  .  0 .  . 

re TUFA 

END 
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In  accordance  with  letter  from  DALN-RDC,  HAHN -AS  I  dated 
22  July  1977,  Suhiect:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Pub l icat ions ,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Curtis,  John  0. 

A  one-dimensional  plane  wave  propagation  code  for 
layered  rate-dependent  hysteretie  materials  /  by  John 
0.  Curtis  (Structures  Laboratory,  U.S.  Army  Engineer 
Waterways  Experiment  Station).  —  Vicksburg,  Miss.  :  The 
Station  ;  Springfield,  Va.  ;  available  from  NTIS,  1981. 

66,  2h  p.  :  ill.  ;  27  cm.  —  (Miscellaneous  paper  /  U.S. 
Army  Engineer  Waterways  Experiment  Station  ;  SL-81-25) 
Cover  title. 

"September  1981." 

"Prepared  for  Assistant  Secretary  of  the  Army  (R&D), 
Department  of  the  Army  under  Project  hAl6ll01A91D." 

Final  report. 

Bibliography:  p.  66 

1.  Electromagnetic  waves.  2.  Finite  element  method. 

3.  0NED3P  (Computer  program),  h.  Soil  stabilization. 

X.  United  States.  Dept,  of  the  Army.  Assistant  Secretary 
of  the  Army  (R&D).  II.  U.S.  Army  Engineer  Waterways 


Curtis,  John  0. 

A  one-dimensional  plane  wave  propagation  code  :  ...  1981 

(Card  2) 

Experiment  Station.  Structures  Laboratory.  III.  Title 
IV.  Series:  Miscellaneous  paper  (U.S.  Army  Engineer 
Waterways  Experiment  Station)  ;  SL-81-25. 

TA7.W3i*m  no.  SL-81-25 
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